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Abstract 

In recent years, bright soliton-like structures composed of gaseous Bose-Einstein conden- 
sates have been generated at ultracold temperature. The experimental capacity to precisely 
engineer the nonlinearity and potential landscape experienced by these solitary waves offers an 
attractive platform for fundamental study of solitonic structures. The presence of three spatial 
dimensions and trapping implies that these are strictly distinct objects to the true soliton solu- 
tions. Working within the zero-temperature mean-field description, we explore the solutions 
and stability of bright solitary waves, as well as their interactions. Emphasis is placed on eluci- 
dating their similarities and differences to the true bright soliton. The rich behaviour introduced 
in the bright solitary waves includes the collapse instability and symmetry-breaking collisions. 
We review the experimental formation and observation of bright solitary matter waves to date, 
and compare to theoretical predictions. Finally we discuss the current state-of-the-art of this 
area, including beyond-mean-field descriptions, exotic bright solitary waves, and proposals to 
exploit bright solitary waves in interferometry and as surface probes. 

1 Introduction 

1.1 Gaseous Bose-Einstein condensates 

In 1925 Einstein predicted that an ideal and uniform gas of bosons, under conditions of sufficiently 
high density and/or low temperature, would begin to "condense" into the single particle quantum 
state of zero energy |TJ. This phenomenon of Bose-Einstein condensation is now known to extend 
beyond gases into liquids and solids, being the underlying mechanism responsible for superfluidity 
in Helium and superconductivity [2]. But it is the gaseous form of this phenomenon that offers the 
purest and most controllable realization of this state of matter []3]-[5j . 

Since 1995, gaseous atomic Bose-Einstein condensates (BECs) have been generated in labora- 
tories world-wide. These gases are extremely dilute, with typical number densities of 10 18 -10 21 
m~ 3 , and the onset of Bose-Einstein condensation occurs at ultracold temperatures of around 100 
nK (3j. Typically, experiments are run sufficiently far below this critical temperature that practi- 
cally all the atoms enter the Bose-Einstein condensate, and the remaining thermal gas component 
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becomes negligible. Within the Bose-Einstein condensate, the individual de Broglie wavelengths 
of the atoms overlap, forming a single coherent matter wave that extends across the system. From a 
theoretical perspective, this enables the description of the many-body system via a single mean-field 
wave equation. 

Although these gases are dilute, the atom-atom interactions play a significant role and intro- 
duce a nonlinearity into the system. At such low temperatures and densities, the interactions arise 
predominantly via elastic s-wave collisions, which are short-range and introduce a local cubic non- 
linearity into the mean-field wave equation. Furthermore, these interactions are usually repulsive. 

The gases are formed and held within confining potentials produced via magnetic or optical 
fields. These make the condensate finite in size and introduce an inhomogeneity across the system, 
both of which have major implications for the static properties and dynamics of the gases, not least 
for the bright solitary waves considered herein. Furthermore, these configurable traps allow for the 
dimensionality of the system to be engineered to produce "quasi-one-dimensional" and "quasi-two- 
dimensional" systems. 



1.2 Solitons and bright solitary matter waves 

Solitons are non-dispersive waves that arise across nonlinear systems, such as shallow water, plas- 
mas and optical fibres [6,7]. Although solitons are defined formally as mathematical solutions of 
nonlinear wave equations, a physical, "working" definition is that a soliton (7J: 

• Retains its initial shape for all time 

• Is localized 

• Can pass through other solitons and retain its size and shape. 

The mean-field wave equation of a BEC is of the form of the (3+l)D cubic nonlinear Schrodinger 
equation, with an additional inhomogeneous term arising from the trapping potential. In the the- 
oretical limit of ID and in the absence of trapping in the remaining direction, this reduces to the 
ID nonlinear Schrodinger equation, for which bright and dark soliton solutions are known to exist. 
Bright solitons are localized humps in the field amplitude, bound together by a focussing nonlinear- 
ity. Dark solitons appear as localized reductions in an otherwise uniform field amplitude, preserved 
instead by a defocussing nonlinearity. While analogs of dark solitons have been observed in BECs 
(see [ 8 ] for a review), we here focus on the case of bright solitons. 

Bright soliton-like (9 - 11 1 structures have been observed in BECs, with the required focussing 



nonlinearity arising from the attractive atomic interactions. Of course, the physical reality intro- 
duces three dimensions and trapping potentials/finite-sized systems, and so these are strictly dis- 
tinct objects from the true bright solitons. As such we will henceforth refer to this wider family 
of structures as bright solitary waves. Following the definition of Morgan et al. fT2| we regard a 
solitary wave as a wavepacket that propagates without change of shape. This relaxed definition will 
allow for the inclusion of solutions which feature trapping potentials and three dimensions, as we 
shall see. Of course, one should not assume that a solitary wave will mimic the wider properties of 
the classic bright soliton and it is a key focus of this Chapter to elucidiate this intimate relationship. 

Bose-Einstein condensates are an attractive system in which to study solitonic waves, with some 
key features summarized below: 

• A sophisticated toolbox based on atomic physics allows almost arbitrary shapes of confining 
potentials to be constructed, for example, waveguides to steer the wavepackets, systems of 
reduced dimensionality, and disordered and periodic potential landscapes. 
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• This toolbox also enables the interactions (i.e. the nonlinearity) to be changed effectively 
from infinitely attractive, through zero, to infinitely repulsive via the exploitation of Fesh- 
bach resonances. Moreover, one can employ atoms such as 52 Cr which feature permanent 
magnetic dipole moments; this introduces long-range atom-atom interactions, i.e. nonlocal 
nonlinearity, into the system [ fT3| . 

• The condensate density can be imaged directly with high contrast. While this is most com- 
monly performed via destructive techniques based on optical absorption, non-destructive 
imaging techniques are also possible, e.g. phase-contrast imaging |3j. The phase of the 
condensate can also be mapped out in space and time via interferometric techniques p4| . 

• Bright solitary waves, which typically exist as small BECs, are mesoscopic quantum systems. 
This scale allows interfacing of the robust and well-established mean-field description of 
BECs with more sophisticated models that incorporate thermal and quantum effects Q. 

• The precision and control offered by BECs makes them an attractive system in general for 
application in ultra-precise force detection and quantum information. For these applications, 
bright solitary waves offer further merits through their self-trapped, highly-localized form. 

1.3 The mean-field Gross-Pitaevskii equation 

Our theoretical analysis will be based upon the well-established Gross-Pitaevskii equation, which 
is a wave equation for the classical field of the many -body wavef unction (3j-^J. This equation is a 
valid description for a gaseous BEC providing: 

• The condensate is macroscopically-populated, i.e. N 3> 1, where N is the number of atoms 
in the condensate. 

• The temperature of the gas satisfies T <^T C , where T c is the critical temperature for Bose- 
Einstein condensation, such that approximately all the particles are within the BEC phase. 

• The dominant inter-atomic interactions are two-body, short-range elastic s-wave collisions, 
whose lengthscale is parameterized by the s-wave scattering length a s . 

• The condition of length scales a s <C d, where d the average interparticle distance, is satisfied, 
such that the detailed shape of the inter-atomic potential becomes unimportant and can be 
modelled by a simple contact potential (hard-sphere collisions). 

• The interactions are weak, parameterized by the condition n\a s \ 2 <C 1, such that fluctuations 
out of the single-particle state are negligible. 

Subject to these criteria, the condensate can be parameterized in time and space by a mean- 
field order parameter i//~(r,t), often termed the macroscopic wavefunction. For convenience we take 
y(r,f) to be normalized to unity, i.e., 

||v(r,0| 2 d 3 r=l. (1) 

According to the Madelung transform, \j/(r,t), which is complex, can be related to the atom number 
density n(r,t) and a phase function 6(r,t) via, 

VKM) = y^pW-0(r,O], (2) 

where ,/V is the number of atoms in the condensate (introduced here to account for the normalization 
of y to unity). 
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The field y evolves in space and time according to the Gross-Pi taevskii equation |3^(5j, 



ih 



dt 



2m 



V 2 + V(r) + 
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m 



V(r,t), 



(3) 



where V(r) specifies the external potential acting on the condensate (taken, for simplicity, to be 
time-independent) and m is the atomic mass. 

The time-independent eigenstate solutions of Eq. ([3]) obey the GPE in its stationary form, 
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where jj. is a (real) eigenvalue. The lowest energy solution to this equation represents the mean-field 
ground state of the BEC. 

The Gross-Pitaevskii equation has proven an excellent description of a vast spectrum of static 
and dynamical properties of BECs [3-5]. The present work will be based primarily on this mean- 



field description, although we will briefly discuss beyond-mean-field descriptions in Section 7. 1 



As is most commonly used to confine a BEC, we will assume a trapping potential that is har- 
monic in shape. For simplicity, we will further assume the trap to be cylindrically symmetric. This 
restriction sacrifices only a little generality for significant gains in clarity. We write this potential 

as, 

V(r) = \m \co 2 x 2 + (0 2 (y 2 + z 2 )] , (5) 



-m[co 2 x 2 + C0 2 (y 2 +z 2 )], 

where co x is the trap frequency in the axial (long) direction and co r is the trap frequency in the 
transverse directions. When dealing with such three-dimensional (3D) systems we introduce the 
trap anisotropy A = (O x /(O r , with X < 1 (> 1) corresponding to a prolate (oblate) trap. 
It is useful to parameterise the interaction strength of the condensate via, 



\a s \N 



(6) 



where a r = ^Jhj mco, is the harmonic oscillator length in the radial direction^] 



1.4 Chapter overview 

The bright solitary waves generated experimentally are related, but strictly different, entities to the 
true bright solitons (which apply only in ID and for a uniform, infinite system). It is the focus of 
this Chapter to explore this relationship in detail, highlighting the similarities and differences. In 
essence, we wish to shed light on how "soliton-like" these solitary waves are. We will consider 
how the waves look, how they move and how they interact with each other. The deviation of the 
bright solitary wave from the true bright soliton is a consequence of two factors: the inclusion of 
an inhomogeneous trapping potential and the extension to three-dimensions, and we will consider 
these two factors separately so as to elucidate their independent contributions to the identity of 
bright solitary waves. We first begin in Section [2] by outlining the experimental generation and 
observation of bright solitary matter waves to date. Following this we begin our theoretical analysis 
of bright solitary waves. Sections [3] and |4] explore the static solutions of bright solitary waves. In 
Section [3] this is conducted within an effective ID model of the condensate, and the role of axial 
trapping considered. Then in Section [4] we extend our analysis of the static solutions to 3D, where 
the collapse instability comes into play. In Sections [5] and [6] we turn to the dynamics of bright 

1 Note that in works that focus specifically on fully trapped condensates, k is more commonly defined in terms of a 
geometric average of trap frequencies (e.g. Refs. [11 15 -24]). The radial harmonic oscillator length here is advantageous 
as it allows us to readily consider the case of zero axial trapping (A = 0). 
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solitary waves in ID and 3D, respectively. There we consider the solitary wave dynamics resulting 
from the presence of axial trapping and the interaction with another solitary wave. In Section [7] 
we turn our attention to the state-of-the-art in bright solitary wave research, discussing beyond- 
mean-field descriptions and the current anomalies with mean-field predictions, observations and 
predictions of more exotic bright solitary waves, and proposals for controllable generation of bright 
solitary waves and exploiting them as ultra-precise atom-optical sensors. Finally, in Section [8] we 
draw some general conclusions. 

2 Bright solitary matter wave experiments 

In order to experimentally realise bright solitary matter waves precise control over the s-wave scat- 
tering properties of an atomic sample is of paramount importance. In the following section we 
discuss the application of magnetic Feshbach resonances as a means of establishing this control and 
review bright solitary matter wave experiments to date. 

2.1 Tuning atomic interactions: Feshbach resonances 

The use of magnetically tunable Feshbach resonances [ 25] to control the interaction between atoms 
is now commonplace in many ultracold atomic gas experiments. Feshbach resonances arise when 
a resonant coupling occurs between the collisional open and closed channels of an atomic system, 
as illustrated in Fig. [TJa). For large internuclear distances, the interaction between two atoms 
can be described by the background potential, Vb g . If two free atoms approach, colliding with low 
energy, E, this potential represents the open or entrance channel for the collision. In contrast, closed 
channels (described by V c ) are able to support molecular bound states. A Feshbach resonances 
occurs when the energy of a quasi-bound molecular state in the closed channel, E c , approaches that 
of the open channel. In this instance a strong mixing between the two channels can occur even in 
the presence of only weak coupling. By changing the magnetic field applied this energy difference 
can be tuned if the magnetic moments of the two channels differ thus the scattering properties of 
the atomic sample can be modified. 

These resonances allow the value of the s-wave scattering length, a s , to be changed over many 
orders of magnitude in both the positive and negative domain by simply changing the magnetic 
field, 

a s {B) = a bg (l - j^-) . (7) 

Here a s is the scattering length at the field of interest, B, ci\j g is the background scattering length away 
from the resonance, A is the width of the resonance and Bo is the resonance position. In the case of 
broad resonances, where A > 1 G, there is a smooth variation of the scattering length through zero 
from positive to negative with a slope of da/dB = a s /A. For Bose-Einstein condensation of some 
species (e.g. Rb, 7 Li) this is of particular importance as it allows the creation of stable condensates 
with repulsive interactions (a„ > 0) despite a negative background scattering length away from the 
resonance. As an illustration, Fig. [TJb) shows the Feshbach resonance in the F = 2,mp = —2 state 
of 85 Rb. This broad resonance, of width 10.7 G, at ~155 G gives tuning of the scattering length on 
the order 40ao/G close to the zero crossing. 

In all of the experiments described in the following sections a Feshbach resonance is the key 
atomic tool without which the controlled creation of bright solitary matter waves would not be 
possible. 
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Figure 1: Feshbach resonances: (a) A two channel model of a Feshbach resonance. A resonance occurs 
when two atoms colliding with energy E resonantly couple to a bound state of the closed channel, (b) The 



Feshbach resonance present in the F = 2,nif = —2 state of 85 Rb 



2.2 Collapse of an attractive Bose-Einstein condensate 

A BEC (in three-dimensions) with attractive interactions is inherently unstable to collapse when 
its interaction parameter k = N\a s \/a r exceeds a critical value k c . This leads to the typical notion 
of a critical atom number N c (for fixed a s and a r ) or critical scattering length a c (for fixed N and 
a r ) at which instability becomes induced. The origin of the collapse instability will be outlined 
theoretically in Section |4~Tj The ensuing collapse of the condensate has been dubbed the 'Bosenova' 
in analogy to the astronomical phenomena of stellar explosion. 



The first experimental insights into BECs with attractive interactions were made using 7 Li |26|. 
Here the negative scattering length of a s = (— 27.4±0.8) an, where an = 5.29 x 10 m is the Bohr 
radius, means that the condensate atom number N grows until it reaches N c and the condensate 
collapses. During the collapse the density of the cloud rises thus increasing both the elastic and 
inelastic collision rates. This causes atoms to be ejected from the condensate with high energy in a 
violent explosion. Following this, the condensate begins to reform, fed by the surrounding bath of 
thermal atoms also present in the trap. If observed for an extended period the system exhibits a saw- 
tooth dynamic of growth and collapse (27J until equilibrium is eventually reached. Throughout, the 
maximum condensate number is strictly limited to the critical number for an attractive BEC. It is 
also possible that collapse occurs even with N < N c due to quantum tunnelling effects and thermal 
fluctuations in the cloud leading to instability. 

Further insight into the collapse phenomena came from the group at JILA (Boulder, US) in 



2001 1 15 28 1, carrying out a controlled collapse using a pure 85 Rb condensate. Tuning the scatter- 
ing length from positive to negative using the broad Feshbach resonance illustrated in Figure [T] not 
only enabled the collapse process to be precisely initiated but also allowed the condition k > k c to 
be fulfilled, unlike systems using fixed negative scattering lengths. Along with control of the initial 
condensate number, control over scattering length made the testing of critical number models pos- 
sible, finding the exact scattering length necessary to collapse the cloud, a c . Early work examining 
the point of collapse using slow field ramps confirmed the relationship between critical number and 
scattering length, determining the critical interaction parameter for k c . Later improvements to the 
calibration of the Feshbach resonance, enhancing precision, found k c to be in excellent agreement 



with mean-field models [29|. 



Following this, the JILA experiments were then extended to study the dynamics of the collapse, 
measuring the evolution of the condensate number following a 'sudden' change in the scattering 
length. Measurements of atom number as a function of time showed a sudden yet delayed loss of 
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Time (ms) 

Figure 2: Controlled collapse: The collapse of a stable Bose-Einstein condensate can be triggered by a 
sudden change of the scattering length from positive to negative. After some time at the new scattering 
length, fcdiapse, the condensate begins to collapse and atoms are lost. Eventually the collapse process ceases, 
leaving a stable remnant in the trap containing Nf atoms. [Data from Durham 85 Rb experiment] 



atoms, as shown in Fig. [2j As the interactions are made attractive the condensate begins to shrink in 
size, thus increasing its density. This contraction tends to accelerate with time eventually leading to 
collapse of the condensate. The time for the collapse to begin t c was found to be shorter for larger 
|a s | as the stronger attraction between atoms in the condensate results in a more rapid contraction 
of the cloud. Following the collapse, a stable remnant component was formed in the trap. Notably, 
the number of atoms maintained in this remnant, Nf, was found to depend strongly on N and a s 
and in many cases exceeded N c . This remnant was observed to persist in the trap for more than 1 s, 
oscillating in a highly excited state. 

In addition, a number of more qualitative features were observed about the collapse process in 
the 85 Rb experiment. The first of these features was bursts of atoms with variable energies being 
ejected from the condensate. These bursts would then focus at multiples of T x /2 and T r /2, where 
T X) r = 2n/(O x ^ r is the trap period in the axial (x) and radial (r) dimensions. In all experiments only 
full, never partial, collapse was observed. However, if interrupted (by jumping the scattering length 
away from the collapse point), jets of atoms were also formed. Unlike the bursts, these streams of 
atoms were found to have highly anisotropic velocities and were interpreted as indicating the local 
pinching of the wavefunction during the collapse. 

The collapse process has since been revisited by the group at the Australian National University 
(Canberra, Australia) |30[ . Again using 85 Rb, measurements of the collapse time have been shown 
to be in good agreement with mean-field models describing the process which take into account 
three-body loss mechanisms. 



2.3 Observation of bright solitary matter waves 

The advent of optical trapping led to the realisation of experimental geometries closer to the ideal 
ID limit. This, in combination with control of the atomic scattering length via Feshbach resonances, 
led to the first observations of bright solitary matter waves by groups at Rice University (Houston, 
US) |9j and Ecole Normale Suprieure (Paris, France) 1 10 1 in 2002 using 7 Li. Despite two inherently 
similar experiments, the ENS group succeeded in producing a single solitary wave whereas the Rice 
experiment resulted in trains of multiple solitary waves. 

In order to utilize the Feshbach resonance in the non-magnetic ally trappable F = I, nip = 1 
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state of 7 Li it is necessary to work using an optical dipole trap Oil] . In both experiments initial 
cooling of the atomic sample was carried out in a magnetic trap using the F = 2,mp = 2 state 
before transferring to a dipole trap and flipping the spin state of the atoms to suppress two-body 
loss mechanisms and allow access to the Feshbach resonance. 

In the ENS experiment optical confinement was realised using a red-detuned crossed dipole trap. 
Here condensates of 2 x 10 4 atoms were produced with a s = +39.7<2o in an approximately isotropic 
harmonic trap. After the creation of the BEC the scattering length was tuned close to a s = before 
adiabatically reducing the power in one of the beams, producing a highly elongated cylindrical 
harmonic trap with co x = 2% x 50 Hz and co r = 2n x 710 Hz. The bias field, and hence scattering 
length, was then ramped to its final value before the vertical beam was switched off, releasing the 
cloud into a ID waveguide. In this trap, the atoms experience a slightly expulsive potential due 
to the magnetic coils used to produce the bias field. As a typical example, at B = 520 G, the trap 
frequency along the waveguide can be considered imaginary, around co x = 2\% x 78 Hz. Tuning the 
scattering length to a small negative value, a s = —3.97 ao, resulted in a soliton of 6 x 10 3 atoms 
able to propagate without dispersion for over 1 . 1 mm. 

In contrast to the crossed ENS trap, the Rice experiment used a single red-detuned dipole beam 
to provide radial confinement. Two additional blue detuned beams were applied to cap the ends of 
the trap in the axial direction. After forming a condensate of 3 x 10 5 atoms with a s « 200<3o the 
magnetic field controlling the scattering length was ramped exponentially to the final value and the 
laser end caps switched off thus setting the resulting solitary waves in motion. 

In this experiment multiple solitary waves were observed. The number of these wavepackets, N s , 
was found to be insensitive to the time constant of the exponential magnetic field ramp. However, 
N s increased linearly with At, the time delay between the switch off of the end caps and the time of 
the scattering length change to a s < 0. For the Rice experiment four solitary waves were observed 
for At = with this number increasing to 10 for At = 35 ms. The wavepackets were observed to 
propagate for ~3 s, this being limited by atom loss rather than dispersion effects. 

With many solitary waves confined in a single trap it becomes possible to explore the dynam- 
ical interactions of the wavepackets. Observation of the solitary wave motion showed evidence 
of a short range repulsive interaction between neighbouring wavepackets raising many questions 
regarding their formation and collisional dynamics. A possible explanation for the formation of 
multiple solitary waves was the presence of a modulational instability [32]. Here, phase fluctua- 
tions of the condensate lead to a local increase in density at wavelengths approximately equal to the 
healing length. The attractive nonlinearity leads to the growth of these density fluctuations and the 
emergence of soli tons. 

The spacing between neighbouring solitary waves observed at Rice increased near the centre 
of the oscillation and decreased near the turning points. This result implied a repulsive inter- 
action between solitary waves. This interaction was attributed to the existence of % -phase dif- 
ferences between neighbouring solitary waves, somehow imprinted during their formation. The 
phase-dependence of the solitary wave interaction will be discussed in Sections 5.1 and[6j and the 



origin of the 7r-phase difference in Section 7.1 



It was not until 2006 that solitary waves were again investigated experimentally, this time at 
JILA 1 1 1 1 using the same 85 Rb experiment that has first observed tunable atomic interactions J33"| 
and controlled collapse [ 15 ]. This new work concluded that the stable remnant observed previously 
in the collapse experiments divided into similar solitary wave structures as seen at Rice. Intrigu- 
ingly, these observations persisted despite the fact that the JILA trap remained almost isotropic (with 
radial and axial trap frequencies of 17.3 Hz and 6.8 Hz respectively), far from the highly elongated 
geometries employed at ENS and Rice. The somewhat surprising capacity of bright solitary waves 
to be supported in almost isotropic trap geometries will be discussed in Section [43] 

Unlike the ENS and Rice experiments, the JILA apparatus used a purely magnetic trap. How- 
ever, the method of creating solitary waves by modifying the scattering length can be considered an 
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Figure 3: Solitary wave oscillation in a weak magnetic trap 1 1 1 1: Following the collapse process used in the 
JILA experiment a stable remnant is formed. The variation in the remnant's width with time can be explained 
by the creation of multiple bright solitary matter waves oscillating in the trap, which are visible in the 2D 
plots of atomic density (insets). 



inherently similar process. After producing condensates of up to 15,000 atoms at a s > the mag- 
netic field was adiabatically ramped to decrease the scattering length to a s = 9ao. To initiate the 
collapse, the magnetic field was changed so as to rapidly (0.1 ms) jump the scattering length from 
a positive initial value to a negative final value of a?. Following some time at the final scattering 
length, £ e voive> the atoms were destructively imaged. Investigating the collapse process as a function 
of Of and the initial condensate number it was clear that the number of condensate atoms surviving 
the collapse could greatly exceed N c . Furthermore, the lifetime of the stable remnant could be as 
long as several seconds. As we will review in Section |6.2| this observation is consistent with the 
presence of several repulsively-interacting bright solitary waves. 

Observations of the condensate size in the trap as a function of time suggested a highly excited 
state had been produced during the collapse, with the remnant cloud's width doubling during its 
oscillation in the trap. However, further analysis revealed that, as in the Rice experiment, multiple 
solitary waves were being created which oscillated back and forth along the weak axial direction of 
the trap, shown in Fig. [3] The wavepackets were observed to persist in the trap for ~ 3 s, undergoing 
as many as 40 collisions in this time. This provided additional experimental data to accompany the 
Rice experiments and the growing body of theoretical work on the stability of three-dimensional 
bright solitary waves (which we will review in Sections [4] and [6]). The number of solitary waves 
created in the 85 Rb collapse experiment was found to be controllable, to a degree, depending on af 
and No. As expected, N s increased with \cif\. Importantly, the number of atoms observed in any one 
solitary wave was never found to exceed N c . 

2.4 Current developments 

In order to further explore the results from both previous experiments and theoretical simulations, an 
experiment has been constructed at Durham University (Durham, UK). As in the JILA experiment, 
this apparatus uses 85 Rb in the F = 2, nip = —2 state allowing access to the 10.7 G wide Feshbach 
resonance giving control over the scattering length of order 40 ao/G close to a = 0. However, 
the trapping geometry, a crossed dipole trap and additional waveguide beam to produce a quasi 
ID geometry, allows entirely independent control of the trap frequencies and s-wave scattering 
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Figure 4: Propagation in the waveguide: (a) As a repulsive BEC travels along the waveguide the interactions 
in the condensate cause it to spread out. (b) In contrast, the attractive interactions present in a bright solitary 
matter wave cause the wavepacket to hold together as it propagates, maintaining its shape with time. [Data 
from Durham 85 Rb experiment] 



properties (due to an independent magnetic bias field). 

Here BECs are first created in the crossed dipole trap (at 300 — 400 «o) by careful tuning of 
the atomic scattering properties. Once condensed, the scattering length is ramped close to a = 
before the BEC is loaded into the waveguide by synchronously switching the cross beams off and the 
waveguide beam on. The scattering length is then jumped to a small, negative value (~ —6.5 ao) and 
the BEC is allowed to propagate along the waveguide as shown in Fig. |4| Weak axial confinement 
along the beam is realised with the addition of a magnetic quadrupole gradient. Typically this results 
in trap frequencies of (O x = 2% x 1 Hz and ft),- = 2% x 27 Hz. Using this method a single soliton can 
be produced containing ~3,000 atoms observed to propagate a distance of ~1 mm in 150 ms. 

In addition to experiments aimed at investigating soliton splitting and binary collisions (the 



theory of which will be detailed in Section 7.2.1 1, the Durham experiment has the potential to be 



extended to the study of atom-surface interactions. Contained within the experimental apparatus is 
a super polished Dove prism (surface roughness < 1 A) designed to allow the study of both classical 
and quantum reflection from a surface. The self-stabilizing, localized nature of the wave packets 
means bright solitary matter waves show great potential as surface probes for the study of short- 
range atom-surface interactions in the future. This idea will be explored in more detail in Section 

The Rice group have reported further experimental activity on bright solitary matter waves 
1 34 1 . Here they '"kick" a bright solitary wave towards a thin potential barrier, formed by a near- 
resonant focussed laser beam. The wave-barrier interaction is observed to result in either reflection, 
transmission and splitting of the solitary wave, depending on the kinetic energy of the solitary wave, 
the potential strength and the nonlinearity (s-wave interaction strength). Moreover, for the case of 
a split solitary wave, they have applied a phase imprinting to one of the solitary waves and thereby 



studied phase-dependent interaction of solitary waves [35]. 



3 Bright solitary waves in ID: static properties 

Having reviewed the experimental formation and observation of bright solitary matter waves to 
date, we will now review our theoretical understanding of these wavepackets. Within a quasi- 
one-dimensional (quasi- ID) system, bright solitary matter waves become completely analogous — 
within the mean-field, Gross-Pitaevskii equation (GPE) treatment — to classical bright solitons of 
the ID nonlinear Schrodinger equation (NLSE) (36j. In this section we examine the quasi-lD limit 
in which this occurs. In Section 13.11 we describe the conventional factorization to reduce the 3D 
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GPE to an effective ID form, and some approaches to include higher-order terms. In Section 3.2 
we demonstrate the link to bright solitons and explore the static properties of bright soliton solutions 
of the NLSE. Then, in Section [33] we consider the bright solitary matter waves which occur as the 
ground state of an axially trapped BEC, elucidating how their form depends on the strength of the 
axial trap, and how they compare to the bright solitons of the ID NLSE. 

3.1 Effective one-dimensional descriptions 

We begin by considering an attractively-interacting (s-wave scattering length a s < 0) three-dimensional 
BEC confined by the cylindrically-symmetric harmonic trap of Eq. ([5]) and described by the 3D 
Gross-Pitaevskii equation, Eq. Q. 

3.1.1 Quasi-one-dimensional GPE 

The quasi-lD limit is associated with highly elongated (co r S> (O x ) traps. The reduction from the 
full 3D to an effective ID description typically proceeds by assuming that the radial confinement 
is sufficiently strong]^] that the radial modes of the condensate become essentially "frozen" into the 
ground harmonic oscillator ground state (i.e. a Gaussian wavefunction) This approximation then 
allows the factorization, 



mco r 

V( r ) = \/^-exp 



-mco r (y 2 +z 2 

2h 



¥ (x) (8) 



where it is implied that both the Gaussian radial wavefunction and the axial wavefunction are both 
normalized to unity. Integrating over the y- and z-directions then yields the quasi- ID GPE for y(x), 



H 2 d 2 m(0 2 x 2 ^ . M , ..; 
-^ 2+ ^--2nco r N\a s \\ V (x, t )\< 



dt 

In the static case, one obtains the stationary quasi- ID GPE, 

h 2 d 2 m(0 2 x 2 | M , . l2 

■^ + — 2 2Hco r N\a s \\ V (x)\ 2 -» 



y(x,t). (9) 



y/(x) = 0. (10) 



This factorization has often been applied in the study of attractively-interacting condensates (in both 
dynamic and static situations) [3"7]442~[ . However, the regime in which this factorization is valid is 
significantly restricted for attractively-interacting condensates [43]; this issue is revisited using a 
full 3D analysis in Section|4] 

3.1.2 One-dimensional equations with 3D effects 

Alternatives to the factorization presented above exist, which yield ID equations retaining more 
3D character by choosing to incorporate the coupling between axial and radial modes, and time- 
dependent dynamics of the radial modes (4^447j. These effects are manifest in the effective ID 
equation through the appearance of higher-order terms. Consequently, the resulting equations have 
a wider range of validity than the bare ID GPE (|9]>, but are no longer isomorphous to the NLSE (for 
co* = 0). 

For example, Salasnich et al. |44,45] chose to factorize the 3D GPE wavefunction into a slowly- 



varying axial function, multiplied by a rapidly varying radial function. The radial function was also 
given a dependence on the axial function itself; this incorporates the effect unique to attractive 



2 Specifically, the criteria fta> r ^> /i and hco r 3> k^T are required to ensure that the condensate and thermal energy 
scales are insufficient to excite the radial modes. 
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interactions in a cigar-shaped trap, where an increase in axial density leads to an associated increase 
in radial density. A variational calculation then yields the non-polynomial Schrodinger equation 

ED. 



in r 

dt 



h 2 d 2 \jf(x,t) 
dx^~ 



+ ■ 



2m 

frCQr 



+ 



m(0 2 x 2 



Y(x,t) + 



2 \ ^l-2\a s \N\¥(x)\ 2 



+ a/1 



2h 2 1 a s \N\ yf(x, t ) | 2 y(x, t) 
ma r ^J\ - 2\a s \N\Y(x,t)\ 2 

2\a s \N\ yr(x,t) \ 2 ) y(x,t) 



(11) 



When |a. s . |A^| V^( x ) | 2 ^ 1 f° r a H x this reduces first to an effective ID equation with both cubic and 
quintic nonlinearities |47 1, and then to the bare ID GPE ([9]) itself. An even more general approach 
can be taken, incorporating even fewer assumptions about the form of the ground state, but leading 
to a coupled system of effective ID equations [46]. 



3.2 Bright soliton solutions 

Consider the ID GPE ([9]) in the homogeneous regime co x = 0. With the removal of the quadratic 
potential term, this becomes a ID nonlinear Schrodinger equation (NLSE) with a focusing nonlin- 



earity [36]. The ID NLSE is a classical field equation which is integrable, in the sense that solutions 
possess an infinite and complete set of conserved quantities [|6[|36j48]. This is analogous to a dis- 
crete system which possesses as many conserved quantities as it does degrees of freedom [49]. This 
integrability leads to a spectrum of true soliton solutions [6,48]. In the case of the ID NLSE with 



focusing nonlinearity, these bright-soliton solutions were first discovered in Ref . [50 5 1 1 using the 
inverse scattering technique (see Refs. (6j|48j for an overview). 

The classical bright soliton solutions of this equation have been extensively studied in the con- 



text of optical solitons [50-56]. The same equation appears in many other fields, including bio- 



physics, astrophysics and particle physics | |36| , and in the study of deep ocean waves |57|. The 
single-bright-soliton solution of the homogeneous ID GPE is given by, 



V(x,t) 



2Jb, 



x exp 



=sech 



a(x — xq —vt) 



2b x 



(12) 



v 2 t 



vix—xo) + — + 



co 2 \a. 



2 N 2 a 2 t 



This solution describes a single bright soliton with amplitude and norrrj^a, velocity v, displacement 
xq, and phase <J>. The parameter b x = %/2m(O r \a s \N is a length scale characterizing the soliton's spa- 
tial extent. Dynamical solutions composed of multiple bright solitons also exist; in these solutions 



each soliton has a similar form to Eq. (12) when well-separated from the others. These multiple- 
soliton solutions contain additional, dynamic phase and position shifts to account for the nonlinear 
interactions between solitons; these dynamical solutions are discussed further in Section|5] 

The single- soliton ground state of the static ID GPE ( 10 1 in its homogeneous (co x = 0) form is 
given exactly by Eq. ( 12 1 with a = 1, v = 0, and arbitrary <I> and xq- The quantity <1> can be chosen 
arbitrarily because it corresponds to a global phase of the wavefunction, and Eq. ( [TO] ) possesses a 
f/(l ) global phase symmetry. Similarly, the displacement xq may be chosen arbitrarily because the 
assumption of homogeneity (co x = 0) ensures the ID GPE to be translationally symmetric. However, 



the choice of displacement xq in Eq. ( 12 1 for the ground state breaks this symmetry; in the context 
of atomic BECs, this symmetry-breaking is a feature of the mean-field description. This feature 
is at odds with a fully quantum-mechanical treatment; in the latter, the ground state of the system 



retains the translational symmetry of the equation, leading to a delocalized ground state [58 59 1. 



3 In contrast to our definition here, a common convention in the literature is to define an amplitude A such that the 
norm is 2A |53| 
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3.3 Effect of axial trapping 



The addition of an axial harmonic trap (co x > 0) removes the integrability of the system and prevents 
the appearance of true solitons. While the new ground state is no longer a soliton, it remains a 
solitary wave in the sense of being capable of propagation without dispersion (see Section [5] and 
Refs. p2][4T][42| ). In this section we elucidate the form of the ground state under axial trapping and 
compare the form of this ground state to the NLSE bright soliton. 



3.3.1 Variational analysis 



A great deal of insight into the form of the ground state can be gained using a variational approach 
p3| ; such approaches have proved useful for treating a variety of problems involving bright solitary 
matter waves | [44| - |47|[60] - |65] ] and will be used extensively in Section [4] 

The ground state solution of Eq. (10) can be alternatively defined as the function \\f{x) which 
minimizes the value of the classical field Hamiltonian, 





~h 2 




j dx 


2m 


lx v{x) 



+ ■ 



m(0 2 x 2 



| \\f(x) | 2 - ha> r N\ a s 1 1 \jf(x) 



(13) 



This functional represents the energy per particle, and generates the ID GPE Q through the func- 
tional derivative SHi D [\j/]/ Sy* = id\j//dt. 

In the homogeneous limit (co x = 0) the ground state is given by Eq. ( [12] ) with a = 1 and v = 0. 
In the trap-dominated limit (co x — > °°) the ground state tends to the harmonic oscillator eigenstate 
y(x) = {m(D x /lih) l l A e~ mmxX ~l 2f '. These limits motivate our use of a normalized Gaussian ansatz, 



y(x) 



m(O x 



1/4 



m(O x x 



exp 



2M 2 



(14) 



x,G 



or a normalized sech ansatz, 



y(x) 



2^/bJ^A 



:sech 



2b x £ 



(15) 



to describe the intermediate regime co x > 0. Substituting theses ansatz into Eq. ( 13 1 we obtain an 
energy functional in terms of the dimensionless length £ x . Note that the sech ansatz length £ Xj $ 
is defined so that £ Xt $ — ^ 1 as (O x — > 0, while the Gaussian ansatz length £ x c is defined such that 
4,g ->lasffl,->oo. 

In the Gaussian case, one obtains the energy functional, 



H lD (£ 



x,G) 




a x \a s \N 
f2%a 2 r £ x £ l 



(16) 



where a x = \Jfij m(O x and a,- = ^Jfij mco r are the axial and radial harmonic oscillator lengths. In the 
sech case, one instead obtains, 



#id(4,s) 



ma> 2 a 2 N 2 



1 



64s 



3f 



1 gjgs ' 
T 24^4 



(17) 



.v.S 



Either of these energy functionals can be analytically (or numerically) minimized to give the corre- 
sponding, variational-energy-minimizing, axial length £ x [43]. The axial lengths £ x for both varia- 
tional solutions are shown in Fig.^a) as a function of the axial trap frequency co x . 



13 




Figure 5: Gaussian-ansatz [Eq. (14i] and sech-ansatz [Eq. ( 15 i] solutions of the ID GPE [Eq. d9b], as 
found in Ref. |43). In (a), the energy-minimizing axial length for each ansatz is shown, as a function of 
(Ox. In (b) the maximum absolute difference between the best-fitting ansatz, (/Ansatz, and the exact numerical 
solution, y/Q, is shown as a percentage of the peak value of \j/Q. This can be expressed mathematically as 
Ay/ = 100max(|y/Ansatz — V / bl)/ max (V / o)- O ur deliberate definition of the ansatz such that £ x $ — > 1 as co x — > 
and £ x ,g — > 1 as co x — > °° results in the potentially confusing trend that 4\G — > as (0 X — > despite the fact 
that the physical length of the Gaussian ansatz tends to a non-zero constant in this limit. 



3.3.2 Comparison to bright soliton solution 

The variational solutions can be compared with a full numerical solution of the ID GPE [Eq. to 
give an idea of how the axial trapping affects the ground state [43 ]. The results of such an analysis 
are shown in Fig. |5jb), which shows the maximum difference in shape between the lowest-energy 
(and hence, most accurate) variational solution and the numerically exact ground state. As one 
would expect, the sech ansatz converges to the exact solution in the axially free limit co x — > and 
it is in this regime, where this ansatz approximates the exact solution well, that the ground state 
can be regarded as soliton-like. In the opposite, trap-dominated limit C0 X — > °°, the Gaussian ansatz 
converges to the exact solution, which is no longer soliton-like in appearance. Convergence is also 
somewhat slower for the Gaussian ansatz as co x — > °° than for the sech ansatz as co x — > due to 
the density-dependent nature of the nonlinearity [43]. In intermediate regimes, one or other ansatz 
provides a good approximation to the solution over a wide range of trap strengths, with only a small 
gap in which neither ansatz is particularly accurate. Consequently, one can usefully think of the 
ground state being deformed from sech-shaped to Gaussian-shaped as co x is increased. 

However, the preceding analysis assumes the validity of the quasi- ID approximation. To ob- 
tain a complete picture of the ground state, and its relationship to the NLSE bright soliton, a full 
treatment of the 3D GPE is required. We undertake such a treatment in the next section (Section[4]). 



4 Bright solitary waves in 3D: static properties 

The 3D GPE [Eq. Q] is non-integrable and does not support true bright solitons. Nonetheless, 
bright solitary matter waves can be observed (9 -[TT | which continue to exhibit soliton-like behaviour 
for a wide range of parameters. They are particularly soliton-like in their dynamical properties — 
especially when considering their mutual interactions and collisions. In this section, however, we 
focus on the static regime by considering the solitary wave ground state solutions of the 3D GPE. 
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The existence and form of bright solitary wave ground states in the 3D GPE is intricately linked 
to the instability of attractive condensates to collapse. In this section we review the properties 
of stationary 3D bright solitary matter waves in detail. In Section 4.1 we discuss the collapse 
phenomena. In Section |4.2| we present variational and numerical approaches to the problem. In 
Section 4.3 we review the properties of bright solitary waves as elucidated by the variational and 
numerical methods, and compare the 3D results to the predictions of the ID approach considered in 
Section |3] 



4.1 Collapse and the critical parameter 



An attractively-interacting BEC in 3D is prone to a collapse instability. Indeed, in the absence 
of trapping, the system will undergo collapse. Importantly, the presence of trapping can support 
metastable, non-collapsing states, although these existence of the metastable state depends on the 
atom number, interaction strength and shape and strength of the trapping potential. The collapse 



instability has been investigated experimentally [15 26-28]. Numerous theoretical studies have 



focused on identifying the parameters associated with the onset of collapse in condensates of various 
geometries, using variational )43][44{[6T||62|[64j , perturbative (24j, and numerical p6|[T7|[23 [|43 W^ 
methods. The condensate dynamics during collapse are the subject of continuing theoretical 



study [30 67-70|. 



Recall, we parameterise the interaction strength of the condensate via k = N\a s \/a r . The rele- 
vance of k is that, when it exceeds a critical value k c , the metastable states cease to exist and the 
collapse phenomenon kicks in. The value of k c is dependent on the trap geometry. 



4.2 Variational and numerical approaches to the static solutions 

The parameter regime of metastable solutions of the 3D GPE with a s < is most accurately deter- 
mined by numerically solving the 3D GPE. However, as shown in Section [33] in ID, a variational 
approach can give insightful and accurate results. Hence we begin with this approach using two 
variational ansatz: an ansatz with Gaussian radial and axial profiles, and an ansatz with a Gaussian 
radial profile and a sech axial profile. 



4.2.1 Variational analysis: Gaussian ansatz 



The solution of the 3D GPE under cylindrically symmetric trapping can be approximated by a 
normalized Gaussian ansatz of the form, 



1 



1/2 



exp 



1 

2a 2 



+ 



v x,G 



I 2 



(18) 



where £ x q and £ r fi are, respectively, axial and radial variational length parameters associated with 
the Gaussian ansatz. (Of course, this becomes the exact solution in the noninteracting regime (a s = 
0). ) Such an ansatz has been considered for bright solitary waves in |43 62|64|, and is most 
appropriate in parameter regimes where the strength of the trap potential dominates over the strength 
of interactions in all directions. Substituting this Gaussian ansatz [Eq. ( 18 1] into the classical field 
Hamiltonian for Eq. ([3]), 
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where V(r) = mco 2 (X 2 x 2 + r 2 ) /2, yields 
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Figure 6: Per-particle energy functional, Z/3D, determined using a Gaussian ansatz [Eq. ( 18 1] for a BEC in a 
cylindrically symmetric, harmonic trap. Trap anisotropies shown are: (a) X = 0, (b) X = 1/2, (c) X — 1, (d) 
X =2, (e) X 2 = —4 x 10~ 4 (expulsive axial potential). The top row [sub-label (i)] shows the case k — 0.35, 
for which all the trap geometries are stable to collapse. In this case there is a stable local minimum in 
the variational energy, which corresponds to the (metastable) bright solitary matter wave ground state. The 
bottom row [sub-label (ii)] shows the case k= 1.1, for which all the trap geometries are unstable to collapse. 



This defines an "energy landscape" in terms of the variational lengthscales £ Xj q and £ r ^, in which the 
variational solution corresponds to an energy minimum. Typical energy landscapes for this Gaus- 
sian variational ansatz are shown in Fig. [6] We seek the lengthscales that minimize this variational 
energy. Differentiating with respect to each of the lengthscale variables produces, respectively, two 
coupled conditions for the variational energy-minimizing lengths, 



and 

*-r,c 



4 1=0. (22) 



2k 
v2n£ x , G 

In the case of prolate and oblate trap potentials these equations can be solved via straightforward 
iterative procedures (43|, while for the axially free case an analytic solution can be found [43 6TJ. 



4.2.2 Variational analysis: sech ansatz 

One can take the same variational approach but with a normalized sech ansatz of the form 

(sdbs) sech (^) exp (-^k) ■ <23) 

where £ x $ and £ r ,s are, respectively, axial and radial variational length parameters. Such an ansatz 
has been considered in [43 6T||62| , and is most appropriate in parameter regimes where the strength 
of the radial trap potential dominates over the strength of interactions, but the strength of interac- 
tions dominates over the strength of the axial trap potential. Following the above procedure, the 
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variational energy expression now becomes, 
H 3D [\j/} = hco r 



1 1 x 2 ^ls , & 



6£ 2 + 2£? + 



24 
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x,S ^r,S " " ^S^S 



yielding the two conditions for the energy-minimizing lengthscales, 



*< S+ ^ S 7T 2 
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The sech ansatz yields variational energy landscapes which are qualitatively very similar to those 
yielded by the Gaussian ansatz [61 62 1. 



4.2.3 Numerical approaches 

A variational approach to the stability of bright solitary matter waves in 3D yields considerable 
qualitative insight, particularly with regard to the collapse phenomenon. However, the approach is 
not particularly accurate in its prediction of the critical parameter k c ; the imposition of a certain 
shape on the wavefunction via the variational ansatz causes variational methods to consistently 
over-estimate k c . Consequently, a great deal of work in the field of attractively-interacting BECs 
has focused on accurately identifying k c , for various trap configurations, via numerical solution of 



the 3D GPE. The main approaches to solving the GPE numerically are reviewed in Ref. |71 1. As 
in Section[3j the numerical and variational results can also be compared in order to investigate how 
bright-soliton-like the bright solitary matter wave ground states are in terms of their shape; such a 
comparison is, however, only meaningful in cases which approach the quasi- ID limit (43). 



Studies have investigated traps with spherical [23]66| and cylindrical |16| symmetry, cylin 



drically symmetric waveguides without axial trapping [62) , and the case of a generally asymmet- 
ric trap IJTtJ . Several works also investigated the configurations of specific experiments in de- 
tail |6Tj|72) . The parameter space of bright solitary wave solutions, under cylindrically-symmetric 
trapping, is summarized in Fig. [7] 



4.3 Static solutions in 3D and the role of trapping 

Here we discuss the predicted bright solitary matter wave solutions (in cylindrically symmetric 
traps) according to the variational method and the full numerical solution. The structure of the 
energy surfaces described by the Gaussian [Eq. ( [20] )] is illustrated for a selection of trap geometries 
and interaction strengths in Fig. [6] The collapse instability is manifest as an unbounded decrease of 
as the variational lengths t x and t r tend to zero. In cases where a bright solitary wave ground 
state exists (upper rows in figures) it is stabilized against collapse by an energy barrier (forming a 
local energy minimum in the energy surface); in cases where such an energy barrier is not present 
(lower rows in figures), no bright solitary matter wave ground state exists. The parameter space of 
metastable ground state solutions as predicted by the variational methods is compared to numerical 
solutions of the cylindrically symmetric 3D GPE in Fig. [7] In this plot we use the parameter A 2 
to specify the trap geometry; this is because that, as well as considering the conventional case 
of confining axial potentials (A 2 > 0) we will also consider the case of expulsive axial potentials 
(A 2 < 0). We will discuss the results of these figures below by separately discussing four key 
trapping regimes (specified in terms of A 2 ). 
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Figure 7: Existence and properties of bright solitary matter wave ground states in cylindrically symmetric 
traps as a function of trap geometry, parameterized by A 2 . The presence of metastable states is indicated by 
the various regions, according to the 3D GPE (blue/grey region), the Gaussian ansatz (region bound by the 
green dashed line) and the sech ansatz (region bound by the red dotted line). The uppermost lines represent 
the critical parameter for collapse k c \ under an expulsive trap A 2 < there exists a lower bounding line 
representing the critical parameter for expansion k e . Note the difference in scale on the abscissa either side 
of A 2 = axis. 



Zero axial potential (A 2 = 0) The case of a zero axial potential (which is equivalent, more 
generally, to any constant uniform potential in the axial direction), results in a waveguide-like trap. 
It leads to some algebraic simplification in the variational equations and, in the case of the sech 
ansatz, the variational energy-minimizing lengths l x and £ r and the critical parameter k c = 3 -1 / 4 
can be found analytically p3||6T| . 

More insight into the physical situation can be gleaned from the corresponding variational en- 
ergy surfaces, shown in Fig.JfJa). The energy surface forms a relatively flat "plain" for larger £ x and 
l r , with sharply rising "ridges" occurring when either length becomes small. However, the (nega- 
tive) interaction term in the energy functional leads to a distinct "chute" [61 ] at the meeting point of 
these two ridges (when both £ x and £ r are small). For low k a raised saddle point separates the chute 
from the plain, thus forming the local energy minimum of the metastable solution; as k increases 
this saddle lowers, until at k = k c it disappears and the entire parameter regime of the plain becomes 
unstable. For the sech ansatz, this transition at k c = 1 /3 1//4 « 0.76 |43 61 1. For the Gaussian ansatz 
the critical value is k c rj 0.778 [64]. For comparison, the non-polynomial Schrodinger equation 
(an extended quasi-lD approach) predicts k c = 2/3, through a simpler calculation [45]. The true 
mean-field result, obtained by numerical solution of the 3D GPE, is k c = 0.675 [62 1. 

Within the regime of metastable solutions, the solitary wave lengthscales vary with the inter- 
action strengths. For k = the axial lengthscale is effectively infinite. As k is increased the axial 
lengthscale reduces monotonically, until the point of collapse. The radial lengthscale stays close 
to the radial harmonic oscillator length a r = \ffij mco r throughout. Interestingly, the solution ap- 
proaches being spherical as the collapse point is reached i foTj . 

In regimes where a bright solitary wave ground state does exist, the energy of the saddle point 
relative to that of the local energy minimum on the plain sets an energy scale at which the bright 
solitary wave ground state will be unstable to collapse when excited. Excitations with sufficient 
energy could allow the condensate to overcome the barrier formed by the saddle point and lead to a 
dynamical collapse in which £ x decreases to zero (6T 62 73[. A second channel of instability also 
arises; because the lack of an axial trap results in a finite- valued energy as £ x — > °°, there exists a 
"dispersive channel" in which excitations of the ground state above a certain energy threshold can 
lead to dynamics where £ x increases without bound | 6*T]62| . 

In Ref. p3j the 3D solitary wave ground state in the waveguide-like trap was compared to the 
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NLSE bright soliton. It was found in Ref. [43 ] that, while it is possible to reach a highly soliton-like 
ground state in a waveguide-like trap, it lies in an experimentally challenging regime. Nonetheless, 
as we review in Section [5] and [6] the dynamics of 3D bright solitary waves can be highly soliton-like 
even when their static shape does not closely resemble the NLSE soliton. 



Prolate and isotropic traps (0 < A 2 < 1) For < A 2 < 1 the trap is prolate, i.e. elongated 
in x, while for A = 1 it is spherically symmetric. In such cases the variational solutions must be 
obtained numerically |43 1. 

The energy landscape under these potentials [with examples shown in Fig. |6jb) and (c)] is 
similar to that for the waveguide trap A = in and around the region of the collapse instability. 
Indeed, the critical point for removal of the metastable state depends quite weakly on A , as evident 
from Fig. [7] The only qualitative difference introduced into the variational energy by axial trapping 
arises in the high-£ A limit, where the potential energy of the trap leads to an infinite total energy in 
the limit i x — > °°, eliminating the dispersive channel altogether. 

In Ref. p3| the solitary wave ground state in a prolate trap was also compared to the NLSE 
bright soliton and similarly to above, the achievement of a soliton-like ground state was found to be 
highly experimentally challenging. 



Oblate trap (A 2 > 1) Such a trapping geometry, in which (O x > (O r , is not typical for the study 
of bright solitary matter waves, as in this geometry no clear analogy can be drawn with an integrable 
NLSE. 

Nonetheless, when an oblate trap possesses a metastable ground state it is indeed a solitary wave 
under the definition used by fT2j . These ground states have been studied using the 3D GPE (62) , 
and 2D reductions with 3D effects |74j. The variational energy surface [Fig. [6jd)] is similar to the 
prolate/isotropic case. 



Expulsive axial potential (A 2 < 0) The self-trapped nature of bright solitary matter waves 
means they can withstand being placed in a trap with a weakly expulsive harmonic axial potential 
(A 2 < 0) without dispersing. This was the case in the experiment of Ref. |10|, and considered 
theoretically in [6T][62| . 

The ensuing variational energy surfaces, shown in Fig. [6fe), again (i) permit metastable states 
[Fig. [fje)i] and (ii) fully collapsed scenarios for k > k c . However, the expulsive potential leads to 
a second instability via an "expansive channel" [62]. This corresponds to axial spreading of the 
solutions i x — > oo. In contrast to the dispersive channel — which never completely prevents the 
existence of a metastable ground state, but renders it unstable to (potentially very small) excitations 
— the expansive channel can destabilize the ground state; like the collapse channel's "chute", the 
expansive channel is separated from the ground state by a saddle point, which disappears for suffi- 
ciently low k, or high | A |. This introduces a critical expansion parameter k e , such that one must have 
k e < k < k c in order to observe a metastable ground state. The structure of k c and k e is illustrated 
in Fig. [7J it is immediately apparent that the regime of metastable ground state solutions with an 
expulsive axial potential is severely restricted compared to the other cases. In particular, | A | must be 
relatively close to zero to avoid passing the cusp point (k c = k e ), beyond which metastable solutions 
are no longer found. 



4.3.1 Asymmetric trap potentials 

Removing the restriction to cylindrically symmetric trap geometries which we have enforced up 
to now leads to a considerably enlarged parameter space to explore. The critical parameter in 
such traps has been numerically determined by Gammal et al. fTTJ. The existence and form of 
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the bright solitary wave ground state in anisotropic traps shows no qualitative differences from the 
cylindrically symmetric case. 



5 Bright solitary waves in ID: dynamics 

When analysing the dynamics of bright solitary matter waves, it is naturally of interest to compare 
their dynamics to the well-known and rich dynamics of bright solitons. The natural regime for such 
comparison is the case in which the axial potential is weak compared to the radial trap potential. As 
discussed in Sections [3] and |4] this regime is where stationary bright solitary matter waves bear the 
greatest resemblance to bright solitons, and experiments to date have focused on this regime. 

Under a quasi- ID geometry, the condensate dynamics are described by the ID GPE [Eq. Q]. 
In the case of zero axial trapping this reduces further to the focusing NLSE, admitting exact bright 
solitons. With a weak trapping or expulsive axial potential, as realized in experiments, integrability 
is lost and the dynamics are no longer those of true solitons. Nonetheless, as we illustrate in this 
section, the dynamics remain highly soli ton-like under the assumption that a ID description is 
accurate. We shall later relax this ID assumption in Section [6] 

In this section, we begin by reviewing the dynamics of single and multiple NLSE bright soli- 



tons (Section 5.1 1, and introduce a particle-like model for their motion and interactions. Introducing 



axial trapping, we then explore the dynamics of bright solitary waves in the quasi- ID approxima- 



tion (Section 5.2.1 ); these dynamics are highly soliton-like and can be easily understood using a 



straightforward modification of the particle model. 

5.1 Dynamics and collisions of the classic bright soliton 

In the absence of axial trapping, an attractively-interacting BEC in a quasi- ID trap is described by 



the focusing ID NLSE, which supports bright soliton solutions [50 51 1. The single- and many 



soliton solutions to this equation have been extensively explored in the context of optical solitons 



[50 -55 1 . We review the results pertinent to soliton dynamics in this section. 



5.1.1 Dynamic bright soliton solutions 

Eliminating the axial trapping in the ID GPE [Eq. yields the focusing NLSE, 

dw(x,t) h 2 d z w(x,t) , ,, , m? , 

ih v y J =-- %±-L±-2h(Q r N\a s • ip(jc,f) \ w(x,t). (27) 

at 2m ax L 



Despite its nonlinear nature, the integrability of Eq. ( 27 ) means solutions can be found using the 
inverse scattering method |50| 51 1. In summary, a scattering transform of yf(x,t) yields, at any time 
t, a spectral decomposition of yf(x,t) into solitons and radiation. The radiation part of the spectrum 
is continuous, and has in general a non-trivial time-dependence. However, the soliton part of the 
spectrum is discrete and time-independent, and is completely described by four real quantities for 
each soliton. Consequently, the spectrum of an /V-soliton solution with no radiation component can 
be completely described by AN real quantities, from which the complete solution y(x,t) can be 
recovered using the inverse scattering transform. 



The most general /V-soliton solution to Eq. ( 27 1, containing no radiation, can be written as 1 53 1, 



where, 



N 

\l/(x,t) = Y,¥j(x,t), (28) 

7=1 



n Y r l + y* i 
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Here we have defined the quantities, 



2 2cOr\a s \N ' 



and, 



Yj = exp 



2 2mco 2 a 2 N 2 
~~h 



(30) 



(31) 



in addition to the characteristic soliton length b x = %/ 2m(O r \ a s \N. Each soliton is described by a real 
amplitude a,-, velocity Vj, position offset xj, and phase <I> ; . In the case that the jth soliton is well- 
separated from the other N —I solitons, the linear system defined by Eq. ( 29 1 can be approximately 
solved to give [53], 



Y(x,t) 




Vjt) 



2b x 



+ 9j 



(32) 



Vj(x-Xj) + ^- + 



(0?a 2 N 2 a 2 t 



Here, qj and are time-dependent position- and phase-shifts which appear as a result of collisions 
with the other N —I solitons. They are given by, 



?j + ?Fj = £±log 

Hi 



dj + a k + i{vj-v k )/2(O r \a s \N 
dj — a k + i{vj — Vk)/2(Q r \a s \N 



(33) 



where the sign is positive (negative) when the jth soliton is to the left (right) of the kth [53 1. While 
the jth soliton is well-separated these shifts remain approximately constant, and only change sig- 
nificantly during collisions. 



5.1.2 Bright soliton dynamics and collisions 

The dynamics of a single bright soliton in the NLSE are determined entirely by their nonlinear 
interactions with the remainder of the solution. It is convenient to divide the remainder of the 
solution into soliton and radiation components, and consider the influence of these components on 
the dynamics separately. We review soliton dynamics due to soliton interactions in this section. 
The majority of these dynamics can be understood on the basis of a simple particle model. The 
interaction of solitons with radiation is more mathematically involved | [52)[53][75| and no similarly 
general picture is available. 

In the absence of radiation, the dynamics of multiple bright solitons are dominated by the inter- 
actions and collisions between solitons. One of the defining characteristics of true solitons, associ- 
ated with the integrability of the system, is that they survive mutual collisions entirely unchanged 
in form. The only observable effects of the collision are the asymptotic position and phase shifts 
discussed in Section l5.1.1l 

The main characteristics of soliton interactions can be illustrated by the collisions of two equal- 
amplitude solitons. This is shown, for various relative phases A<I> = 3>i — <I>2, in Fig. [8] As expected, 
the solitons survive such a collision completely unchanged in form. The position shifts qj are clearly 
visible as the deviation of both solitons from their initial linear trajectories. Although the dynamics 
of the collision itself differ with the relative phase A<t>, the position shift qj is unchanged. Note 
that, due to the phase symmetry of the collision, the 0-relative phase case leads to a central density 
anti-node, whereas for 7i-relative phase a density node is preserved at the origin. 

The independence of the position shifts qj from the solitons' relative phase A<I> allows one, in 
principle, to predict their asymptotic trajectories independently of their phase. Disregarding the 
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Figure 8: Bright soliton collisions in the nonlinear Schrtidinger equation, for soli tons with equal amplitude 
and relative phase A<t> = (a), n/2 (b), n (c), and 3n/2 (d). In each case the density profile of the solution 



is superimposed with the soliton trajectories predicted by a particle model [41 42|. This phase-independent 
model fails to describe the dynamics of the collision in detail, but correctly incorporates the asymptotic 
position shift of the solitons. 



phase information in this way leaves each soliton described by a position, velocity, and amplitude. 
One can then treat the solitons as classical particles with an effective mass proportional to their 
amplitude and some appropriate inter-particle potential. This approach was developed for optical 
NLSE solitons | [76] - |79| , using the inter-particle potential, 

V(xj- Xk ) = ^r^ + r^sech 2 { ^^^) , (34) 

where the solitons are treated as classical particles of effective mass T]y = aj/4. This potential 
reproduces the correct asymptotic position shifts provided the velocities and effective masses satisfy 
the condition |rj ; - — tj^| <C |v; — Vjt|/4a) r |a. s |A^ ||42l. The particle model therefore reproduces the 
asymptotic shift exactly for the equal-effective-mass solitons in Fig. [8] For collisions of solitons 
with non-equal effective mass, the asymptotic shift predicted by the particle model approaches 
the correct value for weak soliton interactions (small density or s-wave scattering length) or short 
interaction times (high-velocity collisions). 



5.2 Dynamics and collisions under axial trapping 

Despite the lack of integrability in the ID GPE with an inhomogeneous axial potential, and the 
resulting absence of true solitons, one may still observe solitary waves if the stationary, eigenstate 
solutions of the equation can propagate without changing shape [12]. While they do not satisfy 
the strict mathematical requirements to be solitons (6j|48j, these non-dispersive solitary waves can 
nonetheless, under certain conditions, behave and interact in a soliton-like way. 



5.2.1 Dynamic bright solitary matter wave solutions 

The possibility to observe solitary waves of this type was examined in considerable generality in 
Ref. [12]. In this work the authors considered, in 1, 2 and 3D, how static eigenstate solutions of non- 
linear Schrodinger equations with a general nonlinearity and an arbitrary external potential behaved 
when used as initial conditions in a nonlinear Schrodinger equation with the same nonlinearity and 
a new, possibly time-dependent, external potential. Two conditions were found to be necessary for 
the original eigenstates behave as solitary waves under the influence of the new potential: firstly, 
the nonlinearity must be decoupled from the absolute position, a requirement immediately satisfied 
by the conventional cubic form of the nonlinearity appearing in the GPE. Secondly, the new po- 
tential must differ no more than linearly in any spatial coordinate from the original potential fT2} . 
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We note that this second condition implies that one can use a time-dependent linear potential as a 
way to control bright solitary waves in experiments without causing them to lose their solitary-wave 
character; experimental control techniques such as this are discussed further in Section [7] 

We shall restrict our attention to the case where the eigenstates are the bright solitary wave 
ground states considered in Section 2, [^] and the "new" potential is a harmonic one of identical 
frequency to the original, but with a displaced center; this is equivalent to the case of the original 
potential acting on a displaced bright solitary wave ground state. In this case, the solitary wave has 
the same spatial profile as the ground state, but its centre of mass moves as a classical particle in 
the static harmonic potential. If free from the influence of other solitary waves or other components 
of the solution, it thus undergoes simple harmonic motion like a classical particle |T2j41]42|. This 



feature of the GPE ground state is analogous to the Kohn theorem for the many-body ground state. 
The latter guarantees that the true quantum mechanical ground state of N bosons in a harmonic trap 
can be expressed as a separable tensor product of a single-body wavefunction in the centre of mass 
coordinate with a general (N — l)-body wavefunction in the remaining inter-particle coordinates. 

5.2.2 Bright solitary matter wave collisions under axial trapping 

If the bright solitary wave is not well-separated from other components of the solution, its dynamics 
are influenced by the nonlinear interaction with the other components. As in the case of bright 
solitons, we concentrate on the interactions between solitary wave components only. For bright 
solitons the asymptotic effects of the interactions were entirely described by phase and position 
shifts (Section 5.1.2| ). For bright solitary waves this is no longer strictly true; however, this can 



be considered a satisfactory approximation in the limit that the external potential is approximately 
constant over the region of the collision and provided that the solitary waves are approximately 
bright-soliton-shaped. Making these approximations, one can combine the particle model of soliton 



collisions (Section [5. 1.2[ ) with the behaviour of a particle in a harmonic trap | |41j|42| . This leads 
to a combined particle model for multiple bright solitary waves in a harmonic trap, which is most 
accurate for (a) weak harmonic traps, (b) fast solitary wave collisions, and (c) in-phase solitary 
wave collisions. 

The collisional dynamics of two, identical bright solitary waves according to the ID GPE are 
illustrated in Fig.[9]for relative phases of and %. As anticipated by the particle model, the dynamics 
are dominated by harmonic particle-like motion when the waves are well-separated; however, when 
the waves collide, periodically, at the trap centre, a soliton-like collision results in a position shift. 



There is no overall phase shift between collisions, however |42|. During the collision, we obtain 
qualitatively the same phase-dependent density structure as for the bright soliton collisions. The 
particle model predicts these trajectories well over the short times shown here. However, over 
longer times deviations do build up, arising from the variation of the harmonic axial potential over 
the characteristic length scale of the collision [ 4T]42| . 



The complex dynamics of three or more solitary waves oscillating and colliding in a harmonic 
trap can be effectively predicted using the particle model; interestingly, the model is itself non- 
integrable for three or more solitary waves, leading to chaotic particle-like dynamics |4Tj42|. 



6 Bright solitary waves in 3D: dynamics 

The solitary waves of the ID GPE explored in the previous section are in many respects similar to 
NLSE bright solitons despite the addition of a trap potential. However, in real experiments it is not 
only the addition of trapping which leads to deviation from the NLSE, but also three-dimensional 



4 It is also possible to consider solitary waves having the form of higher-energy nonlinear eigenstates; such eigenstates 
were considered in Ref . 1 80 1 
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Figure 9: Bright solitary wave collisions in the ID GPE with harmonic axial trapping. The initial solitary 
waves, each with N /2 atoms, are ground states of the ttap displaced by ± w 10.35 in x and with zero ini- 
tial velocity and relative phase A<t> = (a), and n (b). In each case the density profile of the solution is 
superimposed with the soliton trajectories predicted by the particle model. 



effects. While the 3D ground state is still a solitary wave (12) , residual 3D effects can lead to 
large deviations from soliton-like behaviour, although there are regimes where highly soliton-like 
dynamics can still be observed. 

Although the absence of an axial trap potential does not lead to exact soliton solutions in 3D, 
we nonetheless begin by considering the axially free case in Section |6.1| We then consider the 



additional effects of an axial trap in Section 6.2 In broad parallel to the previous section, we 



focus on the dynamics and interactions of solitary waves only. However, this distinction is blurred 
due to the 3D effects during collisions, which can lead to non-soliton-like behaviour and eventual 
destruction of solitary waves. This behaviour is fundamentally linked to the collapse instability in 
3D, explored in Section|4j 

6.1 Dynamics and collisions in a waveguide 

In this section we consider a waveguide-like trap, with harmonic radial and zero axial trapping 
potential. In such a trap the solitary wave profile is the ground state, in the parameter regimes that 
are stable against collapse (Section[4]). With uniform axial potential this solitary wave is self-trapped 
in the .^-direction, and the dynamics of multiple such solitary waves is consequently dominated by 
their interactions, as for NLSE bright solitons. 

6.1.1 Stability of solitary wave collisions 

In the absence of analytic solutions for binary solitary wave collisions in a waveguide trap, such 
collisions must be simulated numerically. This can be done from an initial condition composed of 
two copies of the (numerically obtained) ground state, displaced from each other by some distance 
and given some velocity toward each otheij^] For equal-sized solitary waves the resulting collisions 
can be studied within the parameter space of incident velocity v, interaction strength parameter k, 



and relative phase A<J> [72]. 



5 Such a velocity is imparted numerically by applying a spatially varying phase of e ±,mv -» : A Experimentally, this could 
be achieved by applying a linear external potential to each solitary wave for a short time. 
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Figure 10: Stability of solitary wave collisions in an axially homogeneous waveguide potential, (a) Space- 
time plots of atom density during collision between two identical solitary waves, each with k = 0.4 and 
featuring A<t> = 0, for (i) high incident speed and (ii) low incident speed, (b) The same but for A<t> = n. (c) 
Stability diagram of the solitary wave collisions as a function of incident speed v, and interaction parameter 
k, where the solid (dashed) line marks the boundary between stable and unstable collisions for A<t> = (n) . 
The vertical line denoted k c indicates the critical interaction strength for collapse of a single, isolated bright 
solitary wave. The figure and parameters are from |[72|. 



As for the ground state itself, the key parameter determining the stability of collisions of this 
type is the interaction strength parameter k; this must remain below some threshold k c in order to 
avoid a dynamically-induced collapse when the waves meet. However, k c itself is dependent on the 
other collision parameters. In particular, k c is larger for faster collisions, and for collisions with a 



relative phase closer to K. This is evident from the simulated results in Fig. 10 'a) and (b). The 
latter effect is most noticeable for low velocities, with the phase-dependence of k c disappearing 
in the high-velocity limitj^] At low velocity, this phase-dependence can be understood from the 
collision profiles illustrated for the NLSE in Fig. [8j in the case A<t> = % the density profile of the 
collision itself resembles two solitons interacting repulsively [53j and never overlapping, whereas 
in the case AO = the solitons overlap, leading to a strong density peak. While this peak is of no 
consequence in the NLSE or the ID GPE, in the 3D GPE this peak in the atomic density can trigger 
the collapse instability. The full dependence of the collisional stability on k and incident speed v 
is shown in Fig. 10 for the cases of A<I> = and A<t> = %. Note that the phase dependence of the 
collisional stability is also predicted by effective ID equations retaining more 3D character than the 
ID GPE (47). 

The dependence of k c on the incident velocity v can be understood in terms of the relationship 
between the characteristic time for collapse of the condensate, t c , and the characteristic time for the 
collision-interaction to take place, t- mt . In Ref. | |72| it was illustrated that the critical collision veloc- 
ity, below which collapse occurred in numerical simulations of collisions, for the parameters of the 



JILA solitary wave experiment [11], corresponds to a collision-interaction time ?j nt approximately 
equal to the experimentally measured collapse time, t c . Theoretical investigation of the role of the 
two timescales has not proceeded further to date, in part because the GPE has not been generally 
considered an accurate predictor of t c . However, recent results suggesting that the GPE can accu- 
rately predict t c when a three-body loss term is included |30| offer the possibility of further progress 
in this area. 



6 The GPE is based on the assumption of atomic scattering at low energy and momentum and so by "high velocity" 
here we refer to a scale relative to the condensate's natural speed scale of the speed of sound c = \J 47tha s n/m 2 |3J 
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Figure 11: (a) Population transfer AN during bright solitary wave collisions within a homogeneous waveg- 
uide, as a function of the relative phase between the waves A<t>. At high incident speed (stars) the transfer 
is sinusoidal. For reduced speed (circles) the transfer becomes larger and skewed towards A<I> = 0. For low 
speed (squares) the population transfer diverges as A<t> — > 0, due to runaway nonlinear effects and collapse. 
The results are taken from | |72| . (b-d) Wave dynamics for the low speed case with (b) A<t> = 0, (c) n/2 and 

(d) K. 



6.1.2 Population transfer in solitary wave collisions 

Another effect occurring as a result of the 3D nature of the system is that of population transfer be- 
tween bright solitary waves. In both the ID NLSE, and the 3D GPE for a waveguide trap, collisions 
between solitons or solitary waves with relative phases AO = and % have a density profile which 
remains completely symmetric in x after the collision; in this respect the ID NLSE and 3D GPE are 
analogous. The two descriptions, however, lead to very different dynamics for intermediate phases 
< AO < % and % < AO < 2n. In the ID NLSE the density profile, which is initially symmetric in 
x, loses its symmetry during the collision and regains it afterwards. In the 3D GPE for a waveguide 
trap, the initially-symmetric density profile loses its symmetry during the collision, and this loss of 
symmetry leads to population transfer between the two waves: the first solitary wave grows in am- 
plitude and slows down, while the second wave loses amplitude and speeds up. Example dynamics 
are shown in Fig. [TTJb-d). In addition to the 3D GPE, this effect can also be seen in effective-ID 
approaches retaining extra 3D character [47]. 

The amount of population transferred shows interesting dependencies on the relative phase and 
velocity of the solitary waves (Fig. [IT} [72]. For fast collisions the amount of population transfer 



depends sinusoidally on the relative phase, with the maximum transfer occurring at AO = k/2 
and AO = 3n/2, and the magnitude of this transfer decreasing with velocity. At lower velocities, 
however, this dependence becomes heavily skewed, with the maximum transfer occurring closer 
to AO = and the collapse instability occurring in extreme cases. This deviation from sinusoidal 



transfer is a result of strong transient nonlinear effects during the collision [72]. 



6.2 Dynamics and collisions under axial trapping 

We now introduce an axial harmonic trap, in addition to the waveguide-like trap of the previous 
sections. In analogy to transition from the NLSE to the ID GPE, such a system supports solitary 
wave excitations which have the shape of the system ground state and move like classical particles 
in the axial harmonic trap when well-separated from other solitary waves. When these solitary 
waves do interact and collide, two key questions arise: (a) how soliton-like, and (b) how stable 
are these collisions? Both these questions are particularly pertinent to the interpretation of bright 
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Figure 12: Stability of bright solitary wave collisions in a cylindrically symmetric 3D trap, as established in 
Ref. | |8T| . The number of ID-like bright solitary wave collisions, Cid, is shown as a function of the incident 
velocity v, relative phase <t>, and trap anisotropy A. Here, Cid is defined as the number of collisions for 
which the solitary waves return to within 75% of their original peak amplitude and position at their maximum 
distance from the origin. Higher Cid indicates greater stability. The interaction strength parameter k varies 
with anisotropy as k(X) = ^/A/0.08; this ensures that the effective quasi- ID trap frequency CO — X/4k 2 
remains equal to 0.02 (ST). Rather than being released from displaced positions in the trap, the solitary 
waves are launched with a controlled velocity and relative phase using the interference method described in 
Ref. (STJ. The computation takes advantage of the radial symmetry of the problem, using a pseudospectral 
split-step method in 2D cylindrical coordinates. 



solitary wave experiments to date, and to unresolved issues regarding the role of relative phase (as 
will be discussed in Section[7]). We address each question in turn in the following sections. 

6.2.1 Soliton-like dynamics 

Again, the trapping is assumed to be cylindrically symmetric with a shape specified by the trap 
anisotropy A = (O x /(O r . As was discussed in Section [4] the static bright solitary wave solutions in 
such traps have the most soliton-like shape in the low-A limit (provided the collapse threshold is 
not exceeded), and least soliton-like in the opposite case. One naturally expects a similar trend in 
the soliton-like-ness of the solitary wave dynamics, but is faced with the issue of how to quantify 
the soliton-like-ness of the dynamics. 

When considering repeated collisions between two solitary waves at the centre of a harmonic 
trap, as shown in Fig. [9] soliton-likeness can be defined in relation to the characteristic tendency 
for true solitons to emerge from mutual collisions unscathed. In Ref. | [8T| a metric was defined 
as the number of binary collisions for which the solitary waves subsequently reach their turning 
point in the harmonic trap while still having amplitudes and displacements from the origin above 
75% of their original value. Since ID collisions satisfy this criteria almost indefinitely]^] this metric 
is termed the "number of ID collisions", Cid- While more sophisticated definitions could also be 
employed, this simple metric for the soliton-likeness of the dynamics is enough to reveal a rich 
variation in the dynamics of collisions in the parameter space of incident velocity v, trap anisotropy 



A and relative phase A<I> (81J. Fig. 12 illustrates the variation of Cid in the v-A-A<I> parameter space, 



7 Eventually, the effects of the variation in the external trap potential across the collisions could lead to the break-up 
of solitary waves in the ID GPE. However, this does not seem to occur on timescales easily accessible to numerical 
simulation; instead the numerical errors grow faster than the deviation from the soliton-like behaviour. 
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for interaction strength k = a/A/0.08; this choice of k ensures that the effective ID trap strength 



used in Ref. |81 1 is fixed to CO = 0.02, for which value a solitary-wave ground state exists for N 
atoms. Rather than being set in motion by an initial displacement, the solitary waves in Fig .[72] are 
set in motion using an interference protocol providing arbitrary control over the relative phase and 
velocity (see Section[7]and Ref. (8TJ). 

As would be expected from the analysis of collisional stability in the waveguide trap, how 
soliton-like the solitary wave collisions are is strongly dependent on the relative phase AO at low 
velocity, with Cid being significantly higher for AO = 71 than for AO = 0. This phase-dependence 
weakens for faster collisions, which become more soliton-like as the velocity is increased. How- 
ever, the dependence of Cid on X is oscillatory in character. This is quite distinct from the case 
when considering the static solitary wave ground state (Section [4]), where the solitary wave shape 
varied smoothly with X. These oscillations represent the entirely dynamical effect of radial breath- 
ing oscillations of the solitary waves being excited during the collision. Depending on the trap 
anisotropy X the transfer of energy to these radial breathing oscillations can be be either enhanced 
or suppressed, possibly providing an experimental "knob" with which to enhance the stability of 
bright solitary waves [81]. 



6.2.2 Stability of solitary wave collisions 

While the analysis of the previous section gives a comprehensive account of the soliton-likeness of 



collisions for varying velocity, anisotropy and relative phase, the interaction parameter k in Fig. 12 



is fixed (as a function of X). Other theoretical work has explored the ^-dependence of solitary-wave 



collision dynamics, for fixed X [72 82]. It was found, for example, that even in the presence of 



axial trapping, the stability diagrams in v-k space were qualitatively the same as for the axially- 



homogeneous case [Fig. 10 1. The population transfer for AO ^ 0,71 also occurs within an axially- 
trapped system, but with an additional consequence: the repetition of collisions under axial trapping 
leads to the continued growth of asymmetric populations between the two colliding waves, which 
terminates only when collapse instability occurs in one of the waves. As such even a small deviation 
from AO = or % can, over repeated collisions in a trap, leads to significant changes in the long 



term state of the system 1 82 1. 



These studies also performed modelling of the JILA experiment (TTJ using the 3D GPE. Excel- 
lent agreement with the experiment data was achieved assuming the bright solitary waves to have 
a relative phase very close to AO = 71, i.e. locally repulsive interactions. In this manner, the soli- 
tary waves are predicted to be able to survive many collisions without collapse, as were observed 
in the JILA experiment. This is in general agreement with other work predicting that the multiple 
soliton-trains seen in experiment HfTTJ are consistent with neighbouring solitary waves having a 
relative phase close to 71 [ 22 , 3 8 -40 , 72 , 82 ] . The origin, and even the true existence, of these 7i-phase 
differences is an open question, and will be discussed in Section [7] 

7 Hot topics in bright solitary matter waves 

The degree of control over the potential landscape and nonlinearity of BECs offers unique op- 
portunities to study the fundamental properties of solitary waves, as well as nonlinear systems in 
general. Furthermore, the properties of bright solitary matter waves makes them promising can- 
didates for a variety of future applications. Areas of current research towards future applications 
include the development of soliton atom-lasers 1 83 -85 1, the stabilization and manipulation of bright 



solitary matter waves using spatially and temporally varying traps and inter-atomic s-wave scatter- 
ing lengths [ 86 , 87 ] , and manipulation of bright solitary matter waves in periodic potentials 
with the potential for applications in quantum information (90). 
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In this section we focus on four areas of current research: the description of bright solitary 
matter waves beyond the mean-field approximation (including open questions over their relative 
phase), the exploitation of bright solitary matter waves in interferometry, the application of solitary 
waves as surface probes, and more exotic forms of bright solitary matter wave. 



7.1 Beyond-mean-field treatments and relative phase 

The results presented in this chapter have focussed on the mean-field, zero temperature description 
of a Bose-Einstein condensate provided by the Gross-Pitaevskii equation. However, Bose-Einstein 
condensates are in fact many-body quantum mechanical systems at finite temperature. To incorpo- 
rate these general and important additional effects, more sophisticated models must be employed. 
Such effects become particularly relevant in tightly confined geometries, e.g. a quasi- ID system, 
or close to the transition to Bose-Einstein condensation. The most common family of methods to 
describe beyond-mean-field effects (quantum and/or thermal effects) are the various embodiments 
of the stochastic GPE [4]. A more fundamental quantum mechanical approach is to describe the 
many-body dynamics via the multiconfigurational time-dependent Hartree method for interacting 
bosons (9TJ. 

There have been limited beyond-mean-field studies of bright solitary matter waves to date, al- 
though what studies have been performed have raised intriguing differences from the mean-field 
predictions. 

In the two experiments to date that have generated multiple bright solitary waves (9| 11 1, the 



observed wave dynamics (post-formation of the waves) could be well described by the mean-field 
GPE under the assumption that the waves featured phase differences close to 7T [38 72 82] . This pat- 
tern of relative phases can potentially be explained by a process of solitary wave formation through 
modulational instability, followed by collapse of neighbouring solitons with relative phases close to 
zero [22,38j-40j. However, the mean- field GPE, even when supplemented with phenomenological 



three-body loss terms, cannot provide a quantitatively accurate description of solitary wave forma- 
tion out of a condensate collapse |67j 69 [92j. More sophisticated simulations of the condensate 



quantum field during collapse recover the formation of multiple solitary waves but without 7T-phase 
differences fTOj . The origin of the 7i-phase differences (if they truly exist) remains an open question. 



In Ref. 1 70 1 it was also observed that, under ID quantum field simulations, that the solitary 
wave collisions behaved repulsively, i.e. akin to mean-field collision with 71 -phase difference, but 
independent of the initial relative phase. Although these ID results were not supported by the 
corresponding 3D quantum field simulations, uncertainty remains over whether % -phase differences 
do indeed appear in the experimental systems, and indeed whether relative phase as defined by the 
GPE is a well-defined quantity for experimental bright solitary waves. Potentially, the relative phase 
and velocity dependencies predicted by the GPE could be verified in experiment using the controlled 



generation method of Ref. [81 ], providing a test of the mean-field description. 

More recently, Cederbaum et al. used the time-dependent Hartree method |91 1 to describe the 
collisions between two initially-formed bright solitary waves (93). While the GPE assumes that all 
atoms occupy a single quantum state, this method relaxes this condition and allows occupation of 
an arbitrary number of states. It was observed that the freely evolving solitary waves rapidly lost 
their coherence, evolving into fragmented condensates distributed over multiple states. These are 
fundamentally different objects to the original solitary waves, which within the GPE are assumed 
to perpetuate. These final states can be distinguished in experiment via their signatures in the first- 
order correlation functions. 
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7.2 Bright solitary wave interferometry 



Over the last two decades the advent of atom interferometry [94] has led to significant improve- 
ments in metrological precision for real-world measurements of, e.g., rotation [95] and the accel- 
eration due to gravity [96]. The development of atomic BECs has enabled a new form of atom 
interferometry in which a trapped BEC is coherently split and recombined after a period of differ- 
ential evolution. Following a pioneering early experiment [ |97| , many BEC interferometers have 
been constructed based around the principle of a raised, and subsequently lowered, double-well 
potential [98-102]. This scheme allows long interaction times ^ 1 03] and the small spatial scale 
potentially permits accurate measurements of, e.g., the Casimir-Polder potential of a surface [102 |. 
Provided the raising of the barrier is sufficiently fast, the GPE can provide a good description of 
the dynamics, in the sense that nearly all atoms remain in a single mode, which is coherent across 
the barrier | |104[ . However, interactions also cause undesirable phase diffusion during the interac- 
tion time, and for this reason experiments have typically chosen to reduce or eliminate them where 



possible tl01p03p05 |. 

The properties of bright solitary waves offer a novel solution to the problem of inter-atomic 
interactions; one can envisage an analogue to the optical Mach-Zender interferometer in which 
a BEC is split into two coherent, non-dispersive, spatially-localized bright solitary waves, which 
are manipulated and eventually recombined using a time-dependent external potential. In the fol- 
lowing sections we briefly review existing proposals for the necessary coherent beam- splitting of 
solitary waves, using an internal-state interference protocol [ |8T[ (Section 7.2.1| ) and potential bar- 
riers 1 106f 109| (Section 7.2.2 1. In Section 7.3 we outline proposals to use interferometry devices 
based on bright solitary waves for improved sensitivity in the measurement of atom-surface inter- 
actions [ 1 10]. 



7.2.1 Splitting solitary waves using interference methods 

In this section we consider using a magnetic Feshbach resonance to quasi-instantaneously change 
the (negative) s-wave scattering length from an initial value cq, related to the new value a s by 
a® = a 2 a s . If the BEC is initially in the bright solitary wave ground state at scattering length a° s , 
it remains so immediately after the change to scattering length a s . Assuming a ID description and 
negligible axial trapping (co x = 0), the subsequent dynamics are described by the NLSE [Eq. ( [27] )] 
with initial condition, 

Vb(*) = — ^ sech ( -^Ayr ) > ( 35 ) 



2a^b x \2a 2 b x ) ' 

where the original soliton is assumed to be centered on the origin for convenience. Solutions of 



the NLSE for this initial condition are well-known in the context of nonlinear optics [52]: for 
integer a = J, Eq. ( T35) consists of a bound state, or multi-soliton pulse, of J solitons with unequal 
amplitudes ay and zero velocity (vj = 0). For non-integer a = J + j8, where < j8 < 1, it consists 
of J solitons plus radiation [52]. 

Subjecting such a pulse to a sinusoidal density modulation, such that, 
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:sech 



2a 2 b x 



cos 



/ Kx A<t>\ 

\2a 2 b x + ^2J 



(36) 



is the new initial condition, alters the character of the multi-soliton pulse. For the case of most 



interest, a > 2, this modulation was explored in Refs. |54 55 1 using perturbative and numerical 
methods. In this case, beyond a (relatively low) threshold value of the modulation wavevector K 
the multi-soliton pulse is split into two solitons and a (generally negligible) radiation component. 
The two solitons have equal amplitudes, oppositely directed velocities approximately proportional 



to K [54], and relative phase A<t>. Consequently such a modulation can be used to to control the 
velocity and relative phase of a pair of bright solitons. 
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Figure 13: Interferometric scheme to phase-coherently split a single bright solitary wave into two bright 
solitary waves with controlled velocity and phase |8l) : (a) interferometric protocol illustrated for 85 Rb and 
magnetic field B, which is applied for a short time T between two quasi-instantaneous ^/2-pulses. The 
component in |3, — 2) is not trapped, and escapes, (b) Splitting and re-collision with relative phase (j) —0, and 
(c) = K. Red lines indicate the trajectories predicted by the particle-like model discussed in Section[5]|41 ]. 



In Ref. [8T| , a theoretical scheme was developed to implement such a modulation, and hence 
controllably generate a pair of bright solitons, using an internal state interference protocol; this pro- 



tocol is illustrated schematically in Fig. 13 a). In this protocol, the application of a linear perturbing 
potential for a short time is used to impart momentum on the solitons; the waves continue to behave 
as solitary waves under such a potential as shown in Ref. 1 12]. Numerical simulations verified that 
the same protocol is effective for bright solitary matter waves both in the presence of an axial trap 
[Fig. [T3"]T),c)] and outside the quasi- ID limit. If such a technique can be experimentally imple- 
mented, the possibility to carefully engineer the relative phase between two solitary waves would 
represent a crucial first step towards a general bright solitary wave interferometer. 



7.2.2 Splitting solitary waves at a potential barrier 

Compared to the interference method discussed in the previous section, a simpler method to split 
bright solitary waves is afforded by collisions with a potential barrier. In return for experimental 
simplicity, however, this method lacks the same fine-grained control over the relative phase. 

Within the mean field description, the dynamics of NLSE bright soliton collisions with po- 
tential barriers and wells has been widely explored (see, e.g., | 65| 1 1 If - 114| and Refs. therein). 
The behaviour of solitary waves is similar in soliton-like regimes; in particular, fast solitary wave 
collisions with a narrow barrier lead to smooth splitting of an incoming solitary wave into trans- 
mitted and reflected solitary waves 1 108 109| 114-1 16 1. This behaviour is analogous to bright 
solitons in the NLSE scattering from a 8 -function potential: it can be analytically demonstrated in 
such a situation that the incoming bright soliton is split into transmitted and reflected components, 



each of which consist mainly of a bright soliton, plus a small amount of radiation 1 1 14 1 . Bright 
solitary waves interacting with barriers much narrower than their width largely follow this predic- 
tion 1 108| 115}|116] . Within the mean-field description, potential barrier collisions of this nature 
have been proposed as another means to realize solitary wave interferometers, potentially based on 
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solitary wave molecules | |115| , oscillating bright solitary waves in a harmonic trap (TT6J and bright 
solitary waves in a toroidal trap p 1 08 1. 

In understanding the operation of such an interferometer one must be careful in interpreting 
the predictions of the GPE. This mean-field description is most often thought of as describing a 
system with a wavefunction of Hartree product form (that is, a single macroscopically occupied 
single particle mode). Such a state is free of many-body correlations |107[ . At face value, this 
seems to be at odds with many-body descriptions of the scattering of a bright solitary wave on 
a potential: it was demonstrated in Ref. |106 1, using an effective potential approximation, that a 
condensate bright soliton of 100 atoms could be placed in a coherent macroscopic superposition 
between reflected and transmitted solitons via a slow collision with a wide Gaussian barrier — a 
state entirely dominated by many-body correlations rather than free of them. Similar collisions were 
investigated in Ref. | |107[ using the MCTDHB many-body computational method |91] 117|. In this 
work, the condensate was found to fragment, leaving two macroscopically occupied orbitals. One 
of these orbitals corresponded to a transmitted, and the other to a reflected, bright soliton, implying 
creation of a macroscopic coherent superposition between spatially distinct states. 

This apparent disagreement with the GPE prediction is, however, only a disagreement with the 
Hartree-product interpretation of the GPE; in classical-field methods one commonly uses the GPE 
to describe all macroscopically occupied modes of a system [118], and in this interpretation there 
is no general disagreement with the many-body description. Importantly, in the proposed solitary 



wave interferometers [108 116] the actual interferometric measurement consists of measuring the 
average fraction of atoms ending up on a particular side of a potential barrier. Such a measurement 
is independent of the underlying occupations of single particle modes, suggesting the mean-field 
GPE may still give a useful description of a bright solitary wave interferometer. 

Nonetheless, the realization of macroscopic quantum effects using bright solitary waves offers 
exciting potential for future interferometric devices, as states with macroscopic quantum superposi- 



tion could be exploited to achieve quantum enhancement of the measurement precision [ 103 105]. 
A bright solitary wave interferometer therefore offers the intriguing possibility of observing the 



effects of macroscopic quantum superposition [119 120] through the formation of a fragmented 
state [ 106 107] and hence enhancing measurement sensitivities 1 1 19 121 1. 



7.3 Soliton surface probes and quantum reflection 

There is growing interest in the use of ultracold atoms to measure short range forces close to a 
surface motivated by the possibility of probing short range corrections to gravity which extend 



beyond the Standard model [ 122- 124[ . Traditional experiments, following the pioneering measure- 
ments of Cavendish in 1798 [ |125| , now use a variety of approaches from superconducting gravity 
gradiometers [ |126[ and microcantilevers |127[ to planar oscillators |128[ and torsion balance ex- 
periments pi29] . However, in scaling down experiments to probe ever decreasing length scales a 
new, fundamental problem arises. Quantum electrodynamics predicts a macroscopic force between 
conductors, known as the Casimir force |130| |. This force vastly overwhelms the much weaker 
gravitational attraction between the test masses, such that experiments are forced to search for devi- 
ations between the theoretical and experimental Casimir forces. Precisely calculating such Casimir 



forces for a specific macroscopic test mass near a surface is generally difficult 1 131 1. In contrast, the 



interaction between a single neutral atom and a plane surface is well understood 1 132| 133 1 being 
characterised by the attractive Casimir-Polder potential, 



t/ v dw 

^ret 



c 3 



for z < Aopt/271, 



C4 

--j for A pt/27r < z < A T , 



(37) 
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where for longer length scales the 1/z 3 form of the van der Waals potential, characterised by C3, 
becomes 1/z 4 due to retardation effects. This new regime is characterised by C4 with the transition 
point between the two regimes determined by the wavelength corresponding to the dominant exci- 
tation energy of the interacting atoms, Ao pt 1 134 1. Further from the surface (larger than the thermal 
wavelength of photons, At) the interaction becomes dominated by the thermal fluctuation of the 
electromagnetic field [135]. 

The inherent advantage of directly probing the atom- surface interaction has prompted the recent 
proposal of a new generation of experiments which aim to exploit the precision and control offered 
by atomic physics and ultracold quantum gases to push the measurement of short-range forces into 
a new regime [136 - 139| . Indeed a number of proof-of-principle experiments have already utilised 
ultracold atomic gases to explore the short range van der Waals and Casimir-Polder potentials ]140 - 



143 1 . Nevertheless such experiments are in their infancy and considerable refinement is required 



before they become competitive with the classical 'Cavendish style' experiments as a test of short- 
range gravitational forces. 

The attractive Casimir-Polder potential described above can also be investigated through the 
study of quantum reflection. The term quantum reflection refers to the process where a particle 
reflects from a potential without reaching a classical turning point and is a direct consequence 
of the wave nature of the particle. Significant reflection occurs when the local wave vector of 



the particle k = y (k^ — 2mU(z)/h 2 ) changes by more than k over a distance of l/k, where k^ 
is the wave vector of the particle of mass m far from the potential U (z). This requires an abrupt 
variation in the potential U (z), exactly as is found for an atom in the vicinity of a solid surface. The 
demonstration of quantum reflection from solid surfaces is typically performed at grazing incidence 
in order to reduce the wave vector normal to the surface [ |144[|145] l. The advent of ultracold and 
quantum degenerate atomic samples with large de Broglie wavelengths opens up new possibilities 
to study quantum reflection at normal incidence with unprecedented control over the atomic motion. 
Reflection probabilities as high as 20% have been demonstrated for 23 Na condensates incident on a 
solid silicon surface [146]. 

The use of bright matter wave solitons has been proposed to study quantum reflection from a 
solid surface. In Ref. [ 1 10 1 the authors show that the use of solitons presents a number of unique 
advantages resulting from the presence of attractive interactions. Crucially the robust, self-trapped 
and highly localised nature of bright solitons can result in a clean reflection from the surface, with 
very limited disruption to the density profile as compared to condensates with repulsive interactions 
1 146]. Moreover, previous numerical studies of quantum reflection from purely attractive potential 
wells revealed that in certain regimes the whole soliton reflects with very little loss, leading to a 
significant enhancement of the reflection probability as compared to the single particle case |112| . 
The presence of attractive interactions has also been shown to be advantageous in the performance 
of traps for cold atoms based upon quantum reflection [147|. The absence of dispersion as the 
soliton propagates permits the precise control of the velocity normal to the surface and allows much 
lower velocities to be achieved. 



Current experiments in Durham 1 148 149 1 aim to exploit this combination of advantages and 
promise to deliver accurate measurements of the quantum reflection probability. The proposed 



experimental scenario for the study of quantum reflection is depicted in Fig. 14 



Bright matter wave solitons are formed in an optical waveguide from an 85 Rb condensate with 
attractive interactions. For rubidium and a room temperature surface, the lengthscales relating to the 
Casimir-Polder potential are A opt /27r « 0.12/xm and At ~ 7.6/im. The velocity of the soliton to- 
wards the surface can be controlled by manipulating a weak (~ 1 Hz) magnetic potential along the 
axis of the waveguide. The experimental configuration also allows for the addition of a repulsive 
(or attractive) evanescent field in the vicinity of the surface formed by the total internal reflection 
of a blue (or red) detuned laser field within the glass prism. This produces a potential which de- 
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Figure 14: (a) Schematic of the proposed experimental configuration for the study of quantum reflection of 
bright matter wave solitons from a solid surface. The soliton propagates towards the surface in an optical 
waveguide formed by a focussed 1064nm laser beam. Motion of the soliton along the waveguide is con- 
trolled through the manipulation of a weak magnetic potential along the waveguide. An optional repulsive 
evanescent field can be added through the total internal reflection of a 532 nm laser beam in the prism. The 
inset shows a photograph of the super-polished Dove prism mounted in a UHV glass cell in the Durham 
experiment, (b) The total potential (red) experienced by the atoms in the vicinity of the surface is the sum of 
the Casimir-Polder potential (purple) and the evanescent field (blue). 



cays exponentially with distance from the surface; the decay length being determined by the laser 
wavelength, the refractive index of the prism and the angle of incidence of the laser beam. When 
combined with the atom-surface potential, the repulsive evanescent field leads to a repulsive barrier 



of finite height, in close proximity to the surface (see Fig. 14 (b)). Studies of classical reflection 
from such a barrier can be used to probe the atom-surface potential |140[ . Moreover, the addition 
of both repulsive and attractive evanescent fields can be used to engineer a potential in the vicinity 
of the surface that significantly enhances the quantum reflection probability |150[ . 

7.4 Exotic bright solitary waves 

We have focussed in this chapter on bright solitary waves of a single-species condensate with s- 
wave interactions. However, part of the beauty of these atomic gases is that it is possible to precisely 
introduce additional complexity into the system, e.g. additional condensate components or long- 
range interactions. In certain such scenarios, distinct bright soliton-like structures can arise, with 
even richer properties than their s-wave counterparts. We will briefly summarise these exotic bright 
solitary waves below. 

Bright-dark solitons Condensate systems which involve two co-existing condensates, composed 
of either different atomic species or the same atomic species but in two different magnetic states, 
have been the topic of much experimental and theoretical study [4]. As well as the local intra- 
species interactions within each component, there exists a local inter-species interaction. When all 
interactions are repulsive, and the inter-species interaction is sufficiently strong, phase separation of 
the condensates becomes favourable [3]. Then, the system can enter a state where one component 
adopts a bright soliton-like structure, about which the other component forms a dark soliton-like 
density notch. These hybrid structures, termed dark-bright solitons, have analogs in nonlinear optics 
|151[ and were first predicted to exist in inhomogeneous BECs by Busch and Anglin 1 152 1. Since 



then they have been generated and observed in several experiments [153-155]. Possessing only 



repulsive interactions, they are free from the collapse instability and behave more akin to dark 
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solitary waves (see reference [8] for a review) than bright solitary waves. 



Bose-Fermi solitons Similar to above, it is possible to create a system in which a Bose -Einstein 
condensate co-exists with a degenerate gas of fermions. For identical fermions, the Pauli exclusion 
principle prevents s-wave interactions, and so the predominant interactions in this ultracold mixture 
are the s-wave boson-boson interactions and the boson-fermion interactions. Here, an attractive 
boson-fermion interaction can support a soliton structure in which localized wavepacket of each 
gas, overlapping in space, is self-trapped by their mutual interaction. These structures, termed Bose- 



Fermi solitons, have been predicted via numerical and variational approaches 1 156 - 160], and have 



been simulated to propagate without dispersion [|156[[161]. The boson-fermion interaction must be 



sufficiently strong to overcome the internal repulsion within the BEC component but not so large 
as to induce collapse. As such, in 3D, they exist as metastable states of the system. Interestingly, 
under the inclusion of higher order p-wave interactions with repulsive sign, the collapse instability 
can be completely eradicated, suggesting a greatly enhanced stability |160[ . Related bright solitary 
wave structures are also predicted to arise in Bose-Fermi mixtures in the presence of optical lattice 
potential (162}{T64|. 



Non-local bright solitons In recent years, BECs have been produced in which the atomic species 
have a large natural magnetic dipole moment |13j. In contrast to s- wave interactions which are 
short-range and isotropic, these interactions are long-range and anisotropic. By polarizing the 
dipoles in a common direction, the whole condensate takes on a dipolar nature, which can be ac- 
counted for within the GPE by the inclusion of a non-local term fT3| . The attractive component of 
the dipolar interactions lends itself to support self-trapped solitary wave states, but it also makes the 
collapse instability a general feature of dipolar BECs. 

Within a quasi- ID waveguide, solitary waves of dipolar BECs are predicted to be supported 
[165] 166]. In general, dipolar interactions also co-exist with s-wave interactions, and the capacity 



for self-trapping was shown to occur for various regimes of the dipolar interaction and the s-wave 
interaction. Collisions between two such dipolar solitary waves were found to exhibit more com- 



plicated dynamics including regimes where the colliding waves form a bound state [165] and sen 



sitivity to the polarization direction [166]. More strikingly, in Refs. 1 167[|168[ it was shown that a 



dipolar BEC within a quasi-two-dimensional trap can form a bright solitary wave that is self-trapped 
in two-dimensions and free to move within the untrapped plane of the system. 

Distinct non-local effects can be introduced by coupling of a ground state BEC to highly excited 
Rydberg states ]169) , excited atomic states with high principal quantum number. These excited 
states become coherently shared throughout the condensate, inducing a strong collective interaction 
via van der Waals forces. Using a GPE which incorporated an appropriate non-local term, it was 



predicted in [170| that bright solitary waves could form which are not only free from the collapse 
instability, but are self-trapped in all three dimensions. These wavepackets are predicted to remain 
stable for hundreds of milliseconds and raise the prospect of the first realization of 3D bright solitary 
matter waves. 



8 Conclusions 

We have reviewed bright solitary waves composed of gaseous Bose-Einstein condensates, from 
their experimental formation and observation, through to a theoretical exposition of their static and 
dynamical mean-field properties. Emphasis is placed on how the harmonic trapping potential and 
three-dimensional setting leads to departures in behaviour from the classic bright soliton. Soliton- 
like states remain supported, in the sense that they are capable of self-trapping and retaining their 
static form as they propagate. The most marked deviation is introduced by the extension from 
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ID to 3D, which introduces an instability to collapse for sufficiently large interaction parameter 
and a population transfer during collisions. These deviant behaviours can be greatly reduced in 
appropriate regimes but may be exploited in their own right, e.g. by using population transfer as a 
means to infer relative phase. 

The properties of these bright solitary matter waves are now, generally speaking, well under- 
stood at the mean-field level. Emphasis is now turning to developing a full quantum mechanical 
understanding of these excitations and much work remains to be done in this direction. The experi- 
mental capacity to engineer additional atomic components and interactions into the system promises 
new families of bright solitary waves, which are, for instance, self-trapped in higher dimensions and 
stable to collapse. Furthermore, we believe that bright solitary waves hold strong potential as atomic 
vehicles for applications in matter wave interferometry and surface force detection. 
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